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Abstract 

In this dissertation we study the well-posedness of the three-dimensional La- 
grangian averaged Navier-Stokes (LANS-a) equations. The LANS-a equations are 
a system of PDEs designed to capture the large scale dynamics of the incompressible 
Navier-Stokes equations. In the Lagrangian averaging approach, the motion at spatial 
scales smaller than a chosen parameter a > is filtered without the use of artificial 
viscosity. There are two types of LANS-a equations: the anisotropic version in which 
the fluctuation tensor is a dynamic variable that is coupled with the evolution equa- 
tions for the mean velocity, and the isotropic version in which the covariance tensor 
is assumed to be a constant multiple of the identity matrix. 

We prove the global-in-time existence and uniqueness of weak solutions to the 
isotropic LANS-a equations for the case of no-slip boundary conditions, generalizing 
the known periodic box result [13]. Our proof makes use of a formulation of the equa- 
tions on bounded domains provided by Marsden and Shkoller [21]. In the anisotropic 
model, there are two choices for the divergence-free projection of the viscosity term. 
One choice is the classic Leray projector. In this case, Marsden and Shkoller [22] 
have shown the local-in-time well-posedness of the anisotropic equations in the pe- 
riodic box. We extend their result by considering the second choice of projector, 
the generalized Stokes projector. The local-in-time well-posedness of the anisotropic 
LANS-a with this viscosity term is proven by using quasi-linear PDE-type methods. 

We numerically compute strong solutions to the anisotropic equations in the lam- 
inar channel and pipe by considering steady fluid flow with no-slip boundary con- 
ditions. In particular, given a steady velocity vector that solves the Navier-Stokes 
equations, we numerically calculate the covariance tensor such that the pair solves 
the anisotropic LANS-a equations. Our solutions are in good agreement with the 
results contained in [TTj. Namely, we confirm the logarithmic degeneracy rate of 



the covariance tensor near the boundary and show that in elementary domains the 
sup-norm of the covariance tensor is unbounded in time near the wall. We conclude 
the dissertation by showing the existence of shear flow solutions to the anisotropic 
LANS-a equations. 
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Chapter 1 



Introduction 



1.1 The Navier-Stokes equations. 

The Navier-Stokes equations for an incompressible fluid are a system of partial dif- 
ferential equations that model the velocity vector u = u(t,x) and pressure function 
p = p(t, x) of a fluid whose velocity is divergence- free. We assume the fluid is con- 
tained in a fixed domain Q C lR n with boundary 5Q (possibly empty) and has a 
constant density p with value p — 1. In this case, the Navier-Stokes equations are 
given as 

— — h (u ■ V)u = — gr&dp + uAu + /, 

(J I 

divw = 0, u(0,x) = Uq(x), (1.1.1) 

« = 0on dQ, 
where v > is the kinematic viscosity and the vector / represents external forces. 
The kinematic viscosity is defined as the coefficient of viscosity divided by the density 
of the fluid. The pressure function is determined (modulo a constant) from the 
incompressibility constraint divw = by solving the Neumann problem — Ap = 
div(u ■ V)n with boundary condition gradp • n — 0, where n denotes the outward 
normal vector to the boundary. In the limiting case of viscosity v — > 0, the Navier- 
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Stokes equations are reduce to the idealized setting of the Euler equations. 

The ratio of forces from the convective nonlinearity (u ■ V)w term and the linear 
diffusion uAu term is called the Reynolds number. As the Reynolds number increases, 
the contribution to the motion of the fluid from the viscous term decreases and 
turbulence is introduced. In turbulent regimes, the nonlinear effects send energy 
from the large spatial scales to smaller and smaller scales until the energy reaches 
the Kolmogorov dissipation scale, at which it is abolished by the linear dispersive 
mechanism. To resolve a numerical simulation of the Navier-Stokes equations (jl.l.ip . 
enough grid points or Fourier modes must be used so that the approximation captures 
the energy cascade in all scales down to the Kolmogorov scale. For turbulent flows 
such resolution requirements are not yet achievable, making the problem of turbulence 
an important unsolved problem in physics. The numerical inability in resolving small 
spatial scales motivates the study of the averaged motion of an incompressible fluid. 
The averaged fluid methodology, discussed in the next section, captures the dynamics 
of the large scale motion while averaging the computationally unresolvable scales of 
the Navier-Stokes equations. 

1.2 Averaged fluid motion. 

An approach to modeling the averaged motion of an incompressible fluid is to suppose 
that the velocity of the fluid is a random variable represented by the decomposition 

u(t,x) = U(t,x)+u'(t,x), (1-2.2) 

where U denotes the mean velocity field and vl is a random variable with mean 
zero. In a statistical theory for turbulence, the evolution of the fluid at large spatial 
scales is the primary focus. The process of substituting equation (II. 2. 21) into the 
Navier-Stokes equations and averaging results in the Reynolds averaged Navier-Stokes 
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(RANS) equations [T^ED]- The RANS equations are written as 

dU 

— + (£/• V)U + div(u' (g)u') = - giadp + uAU + f 

L/C 

div U = 0, U(0, x) = u (x) (i- 2 - 3 ) 

U = 0on dn. 
The dynamics of the mean velocity U is well-defined when the Reynolds stress tensor 



(u' <g> w') is expressed in terms of [/. The effect the motion at small length scales 
has on the evolution of the averaged velocity U is called the turbulence closure 
problem. Classically, it is assumed that the Reynolds stress term is of the form 



div [u 1 ® u') = ve^S", x, Def U) ■ Def U, where ve is the eddy viscosity and Def U is the 
rate of deformation tensor defined as 



DeiU = - 
2 



VU+(VUy . (1.2.4) 



Under this assumption, viscosity, that is not naturally present in the physical model, 
is added into the system. This artificial viscosity augments the inherent dissipative 
mechanism and assists in the removal of the energy contained in the small scales at 
which v! resides. Since it is still necessary to guess the form of z/g, an improvement 
to the procedure of modeling the averaged motion of a fluid is needed. 

The Lagrangian averaged Navier-Stokes (LANS-a) equations are a system of par- 
tial differential equations designed to capture the large scale dynamics of the incom- 
pressible Navier-Stokes equations without the use of artificial viscosity or dissipation. 
In the Lagrangian averaging approach, the motion at spatial scales smaller than a 
chosen parameter a > are filtered. The inviscid form of the LANS-a equations, 
the Lagrangian averaged Euler (LAE-a) equations, first appeared in Holm, Marsden, 
and Ratiu [TTI [T6] as a n-dimensional generalization of the one-dimensional Camassa- 
Holm equation. The authors expressed the LAE-a equations as the Euler-Poincare 
equations corresponding to a Lagrangian given by an H ^equivalent norm. In the 
next section, we discuss their results in greater detail. 
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Marsden and Shkoller provide a complete derivation of the LANS-a equation in 
[2~T| [22]. Rather than averaging at the level of the Navier-Stokes equations, the au- 
thors average at the level of the action functional. For the LAE-a equations, the 
action functional is an a 2 -modification to the kinetic energy action functional of the 
classical Euler equations. Unlike the Reynolds averaging procedure described above, 
averaging at the level of the action functional preserves the variational structure of 
the incompressible fluid. In particular, the solutions to the LAE-a equation are the 
extrema of an H ^equivalent energy functional just as the solutions to the Euler equa- 
tions are minimizers of the total kinetic energy. To highlight additional features of the 
Lagrangian averaging approach, we outline Marsden and Shkoller's recent derivation 
[22]. For a chosen parameter a > 0, the authors define the initial data u e = u + eu>, 
w 6 S 12 , in a ball of radius < e < a centered at given initial data uq. Solving 
the Euler equations with initial data Uq and u e Q results in a velocity field u and a 
perturbed velocity field u e . Define r\ and if as the Lagrangian trajectories associated 
with the velocity solutions u and u e . The vector rj solves the first order initial value 
problem 

fj{t,x) =u(t,r)(t,x)), 

(1.2.5) 

77(0, x) = x, 
where the dot represents the partial derivative with respect to time. The vector if 
solves equation (I1.2.5P with u replaced with the solution u e . We can regard the differ- 
ential equation (1 1 . 2 . 5 [) as a map from the space of u's (spatial or Eulerian description) 
to the space of (77,77) (material or Lagrangian description). The Lagrangian fluctua- 
tion £ e := if o r/ _1 is a volume preserving diffeomorphism that plays the role of the 
Reynolds decomposition (II .2.21) . Since the decomposition is made in the Lagrangian 
reference frame, the interplay between the Euler and Lagrangian frameworks is of 
fundamental importance. This interplay is insignificant in the Reynolds averaging 
approach where the decomposition and the averaging only occur in the Eulerian ref- 
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erence frame. After asymptotically expanding u e in terms of e, Marsden and Shkoller 
take the ensemble average of the total kinetic energy 

S=- [ [ \u £ \ 2 dxdt 
2 Jo Jn ' 

over all possible solutions u e . At this point, the turbulence closure problem is con- 
fronted. In the Lagrangian averaging approach, the turbulence closure problem 
amounts to specifying the fluctuations £' and £" ( ' denotes the derivative with respect 
to e evaluated at e = 0) as functions of the mean velocity u. The generalized Taylor 
"frozen turbulence" hypothesis provides the solution and assumes the fluctuation £' is 
frozen or Lie advected in the mean flow as a divergence-free vector field. In addition, 
the fluctuation £" satisfies 

where ( • ) is the ensemble average over the solutions u e and the operator -^ is the 
usual total derivative. Marsden and Shkoller's derivation results in an averaged action 
functional which includes all terms up to order a 2 . To conclude their calculation, the 
authors apply Hamilton's principle to the averaged action functional therefore yielding 
the LAE-a equations. 

The derivation outlined above has recently been generalized by Bhat et. al [I] 
in their computation of the LAE-a equations for a compressible fluid. Unlike the 
previous derivation, the authors averaged over a tube of trajectories centered around 
a given Lagrangian flow. The tube is constructed by specifying the Lagrangian fluc- 
tuation C, £ at t = and deciding on a "flow rule" which evolves £ £ to later time. For 
example, the flow rule for an incompressible fluid is the frozen Taylor hypothesis. In 
the model of compressible fluid motion, the flow rule is the answer to the closure 
problem and is chosen to be one of two different physical properties. One physical 
property is appropriate for isotropic fluid conditions and the other property yields 
the anisotropic model for bounded fluid containers. The remainder of the author's 
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construction follows the outline presented in the previous paragraph. 



1.3 Geometric structure of the Euler and LAE-a equations. 



The Euler-Poincare equations, developed originally by Poincare |26] in his study of 
Euler-type equations, are determined once a Lagrangian map L : q — > R is specified 
in a Lie algebra q. For any point ( G Q, the evolution of the variable ( is determined 
by the Euler-Poincare equations 

d SL . 5L , 

lur^rc (L3 - 6) 

The map ad^ : g — > q, the adjoint representation of the Lie algebra, is the linear 
map t] i— >• [C,rj\, where [C,v] denotes the Lie bracket of ( and t]. The map ad£ 
is the dual linear map associated with ad^. The LAE-a equations and the classical 
Euler equations are the Euler-Poincare equations posed on the same Lie algebra using 
different Lagrangian maps. The configuration space that yields the appropriate Lie 
algebra for incompressible fluid motion was unknown until Arnold [2] and Ebin and 
Marsden [12]. 

Ebin and Marsden [12] define the configuration space for incompressible fluid 
motion as 

X>*(fl) := {ri(t, •) : fi -»• Q\r) E H s (tt),r] is a bijection, 

r]~\t, ■) : Q^n is in H S (Q), anddet Dr) = 1}, 

the group of H s -c\a.ss volume preserving diffeomorphisms on Q. The group V s (Q) 
(under composition) is a C°° differentiable manifold but is not a Lie group [12] (right 
multiplication is smooth, but left multiplication is not). It does however have an 
exponential map and associated Lie algebra in the usual sense of Lie groups. The 
tangent space to T> S (Q) at the identity is identified with the space Xl iv (Q), the space 
of H s divergence-free vector fields on Q that are tangent to the boundary dQ. 
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The Euler equations can be written as the Euler-Poincare equations posed on the 
Lie algebra A^ v (fl) with the Lagrangian defined as the total kinetic energy. This con- 
nection was first made by Arnold [2J. The Euler equations arise from an application 
of Hamilton's principle of least action to the L 2 Lagrangian 

£( u ) = - \u{t,x)\ 2 dx. (1.3.7) 

2 Jq 

Euler-Poincare reduction techniques (see Marsden and Ratiu [20]) show that Hamil- 
ton's principle reduces to the following variational principle with respect to Eulerian 
velocities: 

5 / £{u)dt = 0, (1.3.8) 

J a 

which should hold for all variations 5u of the form 

5u = w + (u ■ V)w — (w ■ V)u, 

where w is a time dependent vector field representing the infinitesimal particle dis- 
placement vanishing at temporal endpoints. The Hamiltonian structure, along with 
references to the literature, can be found in Marsden and Weinstein [23], Arnold and 
Khesin [3], and Marsden and Ratiu [20] . 

Arnold [2] discovered that solutions to the Euler equation correspond to geodesies 
of T> S (Q) with respect a L 2 right invariant metric. The L 2 metric is defined to be the 
weak Riemannian metric on V s (Q) whose value at the identity is 



(u(t,x),w(t,x)) L 2 = / u(t,x) ■ w(t,x) dx, u, w E X^ V (Q). (1.3.9) 

Jo. 

The word weak is used here because this metric need not define the topology on 

the tangent space but may define a weaker topology. The bridge between the group 

V s (Q) and hydrodynamics is the following. If r](t,x) E V s (Q) is a geodesic with 

respect to the L 2 right-invariant metric (II .3. 91) and the velocity of the fluid is defined 

as u(t,x) = f]{t,ri~ l {t 1 x)), then velocity vector u is a solution to the classical Euler 
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equations 



du . _ N , „ 

— - + (it ■ V)u = - grad p + f, 

dt (1.3.10) 

divw = 0, u(0,x) — u(0), 

with the boundary condition that u is tangent to dVL. The geodesic r), which satisfies 
the first order differential equation (11.2.5p . is the particle path associated with the 
vector velocity solution to the Euler equations. By showing the existence of a smooth 
geodesic on the group T> S (Q) for s > (n/2) + 1, Ebin and Marsden [12] obtain short- 
time well-posedness of the Euler equations on a smooth n-dimensional Riemannian 
manifold. 

We summarize the geometric results of the Euler equations on Q with the following 
theorem. 

Theorem 1.3.1 (see |16j ) The following statements are equivalent: 

(i) The velocity vector u £ Xl iv (Q) solves the Euler equations U.3.10\) . 

(ii) The particle path 77 G V s (VL), given by U.2.5\) , is a geodesic on T> S ^(Q) 
with respect to the L 2 right invariant weak metric U.3.9\) . 

(Hi) Hamilton's principle of least action U.3.8\) holds for variations of the 
form 5u = w + (u ■ V)w — (w ■ V)u. 

(iv) The Euler- Poincare equations 

d_M__ M 
dt 5u u 5u 



holds for the Lagrangian £ defined by \1.3.1^ on the Lie algebra X^ V (Q) 



The Lagrangian averaged Euler equations contain a rich geometric structure sim- 
ilar to the geometric framework of the Euler equations. The LAE-a equations were 
first written by Holm, Marsden, and Ratiu [16] as an n- dimensional generalization 
of the one-dimensional Camassa-Holm equation. For a given parameter a > 0, the 
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authors define the Lagrangian for the mean fluid velocity u G Af| iv (Q) as the H 1 - 
equivalent norm 

L(u) = - J u-u + 2a 2 Befu:Beiudx, (1.3.11) 

where the differential operator Def is the deformation tensor defined by (jl.2.4p and 
: is the contraction on two indices given by a : b = dijbij. This H 1 Lagrangian de- 
fines a right-invariant weak metric on the Lie algebra of divergence-free vector spaces 
A*j- V (fi). Holm, Marsden, and Ratiu [16] calculate the Euler-Poincare equations with 
the Lagrangian defined by equation (11.3. 111) . The result is the Lagrangian averaged 
Euler equations in Euclidean space 

f)l) rp 

(it ■ V)w+ [Vm] • v — — grad p, 



dt 

v = (l-a 2 A) Uj (1-3-12) 

div v = 0, v(0, x) = v , 
with the boundary condition that v is tangent to dQ. The authors conclude that the 
solutions to the LAE-a equations correspond to the geodesies of T> S ^(Q) with respect to 
the H 1 right-invariant metric defined by (11.3.111) . In [27] and [28], Shkoller generalizes 
equation (11.3.121) to a smooth n-dimensional Riemannian manifold. By showing the 
existence of smooth geodesic flow on T> S (Q) with respect to the H l metric defined by 
(11.3. lip . Shkoller proves the smooth-in-time local (global when n = 2) existence and 
uniqueness of strong solutions with H s initial data, s > in/2) + 1. 

1.4 Organization of the dissertation. 

In this dissertation we study the well-posedness of the Lagrangian averaged Navier- 
Stokes equation on fluid containers located in M 3 . There are two types of LANS-a 
equations: the anisotropic version in which the (fluctuation) covariance tensor is a 
dynamic variable that is coupled with the evolution equations for the mean velocity, 
and the isotropic version in which the covariance tensor is assumed to be a constant 
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multiple of the identity matrix. A brief history of the LANS-a equations is presented 
in Chapter [2J 

In Chapter [21 we consider the isotropic LANS-a equations on bounded domains. 
We prove the global-in-time existence and uniqueness of weak solutions for the case of 
no-slip boundary data, extending the periodic box result of Foias, Holm, and Titi |14j . 
We make use of a formulation of the LANS-a equations on bounded domains given 
by Shkoller [28] and Marsden and Shkoller [5T] which reveals the correct boundary 
conditions. We begin the proof of existence with a sequence of strong solutions to 
the LANS-a equations. Known to exist from [21], these vector fields trivially solve 
the weak formulation of the equations. We use standard interpolating inequalities 
to establish estimates independent of the sequence index. Classical compactness 
arguments enable us to conclude that the sequence converges in the correct space 
with the limit satisfying the weak formulation of the LANS-a equations. At the end 
of the chapter we prove the existence of a nonempty, compact, convex, and connected 
global attractor. 

Chapters [3] and H] focus on the anisotropic equations. In bounded domains, the 
covariance tensor plays a prominent role in the mechanics of fluid motion. The 
anisotropic LANS-a are therefore the correct model to capture the large scale motion 
of the fluid. There are two choices for the divergence-free projection of the viscos- 
ity term. One choice is the classical L 2 -orthogonal Leray projector. In this case, 
Marsden and Shkoller [22] show that strong solutions exist and are unique in the 
three-dimensional periodic box for a finite time interval. In Chapter [3j we extend 
this result by considering the second choice of projector, the generalized Stokes pro- 
jector. The generalized Stokes projector, defined in detail in Chapter [3], assigns to 
the divergence-free vector field the no-slip boundary condition. Without fluid motion 
on the boundary, the fluctuation tensor is also zero on the boundary, making the 
Stokes projector the appropriate projector for the anisotropic equations. The inclu- 
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sion of the generalized Stokes projector viscosity term in the anisotropic equations on 
a periodic box is an essential step towards understanding averaged flow on bounded 
domains. We prove the local-in-time existence and uniqueness of classical solutions to 
the anisotropic LANS-a on the periodic box with this viscosity term by using quasi- 
linear partial differential equation type methods. We begin the proof by obtaining 
an approximate solution using the Galerkin projection of the anisotropic equations 
onto a finite dimensional vector space. The generalized Stokes projector term forces 
us to repose the problem in terms of the momentum rather than the velocity. After 
proving an elliptic regularity- type result, we show that our approximations remain in 
the correct vector spaces independent of the projection. We use classical compactness 
arguments to establish the existence of solutions for a short period of time. 

In Chapter HJ we examine the anisotropic LANS-a equations in the channel and 
pipe geometry. Recently, Coutand and Shkoller [TT] have proposed a turbulent chan- 
nel flow theory founded on the anisotropic equations. By posing the problem in the 
correct functional framework, Coutand and Shkoller [TT] proved that weak solutions of 
the anisotropic LANS-a equations exist and are unique throughout the entire channel 
for all time. In Chapter m we numerically compute strong solutions to the anisotropic 
equations by considering steady fluid flow with no-slip boundary conditions. In par- 
ticular, given a steady velocity vector that solves the Navier-Stokes equations, we 
numerically calculate the covariance tensor such that the pair solves the anisotropic 
LANS-a equations. Our solutions are in good agreement with analytic results con- 
tained in [TT] . Namely, we confirm the logarithmic degeneracy rate of the covariance 
tensor near the boundary and show that in these elementary domains the sup-norm 
of the covariance tensor is unbounded in time near the wall. We conclude the chapter 
by showing the existence of a shear flow velocity solution to the anisotropic LANS-a 
equations. 

This research was partially supported by the National Science Foundation under 
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Chapter 2 

Weak solutions of the isotropic 
LANS-a equations 



2.1 Introduction. 

The isotropic Lagrangian averaged Navier-Stokes (LANS-a) equations for an incom- 
pressible viscous fluid moving in a bounded fluid container Q C M n with smooth (at 
least C 3 ) boundary dQ may be written as the following system of partial differential 
equations: 

d t u + V u u + U a (u) = -(1 - a 2 A)- x gr&dp - vAu + F , (2.1.1a) 

divn = 0, (2.1.1b) 

u = 0ondn, (2.1.1c) 

u{0,x) = u {x), (2.1. Id) 

where 

U a (u) = a 2 (l - a 2 A)^ 1 Div(Vn • Vu T + Vu-Vu- Vu T ■ Vn), (2.1.2) 
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and d t represents the partial derivative with respect to t. We use u(t,x) to denote 
the large-scale (or averaged) velocity field of the fluid, assumed to have constant 
density. The pressure function p(t, x) is determined (modulo constants) from the 
incompressibility constraint (12.1. lbl) . The constant v denotes the kinematic viscosity 
of the fluid, and a > is the spatial scale at which fluid motion is filtered, i.e. 
spatial scales smaller than a are averaged out. The additional term T{x) G H 2 C\Hq 
represents the external force acting on the system and is assumed, for simplicity, to 
be time- independent. We let A := —PA denote the Stokes operator, with P the 
Leray projector onto divergence-free vector fields. 

It has been a longstanding problem in fluid dynamics to derive a model for the 
large scale motion of a fluid that averages or course-grains the small, computationally 
unresolvable, scales of the Navier-Stokes equations. The LANS-a equations provide 
one such averaged model, and have been studied rather extensively from both the 
analytical, as well numerical, points of view. The numerical simulations of both 
forced and decaying isotropic turbulence given by Chen et al. [51 El El IB] and Mohseni 
et al. [23] are briefly reviewed in Chapter |U 

The LANS-a equations can be written in an isotropic and anisotropic version. 
The isotropic equations, given by equations (I2.1.ip are the focus of this chapter. As 
discussed in Chapter [1], the inviscid [y = 0) version of equations (12.1.11) . known as the 
Lagrangian averaged Euler (LAE-a) or Euler-a equations, was first given by Holm, 
Marsden, and Ratiu [16] in the case that Q = M. n as 

d t v + V u v + [Vwl ■ v — — grad p, 

1 J (2.1.3) 

divw = 0, 
where the variable 

v = (1 -a 2 A)u 

may be thought of as the momentum. Foias, Holm, and Titi [TJjJ first added viscous 
dissipation to (J2.1.3P ; they argued on physical grounds that the momentum v rather 
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than the velocity u, need be diffused. By assuming periodic boundary conditions, 
they obtained the following form of the LANS-a equations: 

d t v + V„f — a 2 \Vu] ■ Au = — gradp + vAv + g, 

divw = 0, (2.1.4) 

u(0,aO = u (x), 

with g taken in L 2 . While yielding the correct equations on a periodic box, the 
question of how to appropriately prescribe boundary data in the no-slip u = case 
remained open. Specifically, inversion of the dissipative term uAv = u(l — a 2 A)u, a 
fourth-order operator, requires further constraints than simply u = 0on dQ. 

Shkoller [28] and Marsden and Shkoller [2T] supplied the additional boundary con- 
dition by reformulating ( 12. 1.4ft as the system of equations ( 12. 1.1ft . In the formulation 
(12.1.11) . it is clear that if u = on dQ, then 

Au = on dQ 

as well. This follows since each term in the inviscid equations identically vanishes on 
the boundary, thanks to the inversion of (1— a 2 A) with Dirichlet boundary conditions. 
The viscous term in the formulation (12.1.11) was obtained by treating the Lagrangian 
trajectory as a stochastic process, and replacing deterministic time derivatives with 
backward-in-time mean stochastic derivatives, exactly following the usual procedure 
for obtaining the viscous dissipation term in the Navier-Stokes equations as done by 
Chorin [9] and Peskin |25j . 

The term U a (u), given in equation (12.1.21) . provides a regularization to the Navier- 
Stokes equations which is dispersive, rather than dissipative, in character. This reg- 
ularization is geometric in nature, and arises as the geodesic flow of an H 1 right- 
invariant Riemannian metric on the Hilbert group of volume-preserving diffeomor- 
phisms of the fluid container (discussed in Chapter [1]); as such, this regularizer yields 
an a priori L°° — H 1 estimate in three-dimensions (in general, n dimensions for n > 2), 
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and one is thus tempted to ask whether the LANS-a system is globally well-posed 
(even though, as is clear from (j2. 1.1 aft , no additional artificial viscosity is being added 
to the Navier-Stokes equations). 

In the case of the periodic box, Foias, Holm, and Titi [2] proved the global 
well-posedness of H 1 weak solutions in dimension three, but as we noted, their for- 
mulation (12.1.41) did not provide the obvious extension to bounded domains. Using 
the equations (12.1.11) . Marsden and Shkoller [2T] proved the global well-posedness of 
classical solutions in dimension three in the case of no-slip boundary data. In this 
chapter, we give an extension of that result to the H 1 weak solutions of Foias, Holm, 
and Titi [2]. In particular, we prove the global-in-time existence, uniqueness, and 
regularity of weak solutions to the LANS-a equations for initial data in the class 
{u G Hq\ divw = 0}. The analogous two-dimensional result follows trivially (as it is 
already known for the original Navier-Stokes equations). 

The well-posedness result leads to the existence of a nonempty, compact, convex, 
and connected global H l attractor in both two- and three-dimensions. In three- 
dimensions, the global attractor has the identical bound as that obtained by Foias, 
Holm, and Titi [2] for periodic boundary conditions. This upper bound depends 
on 1/a and consequently tends to infinity as a tends to zero. Due to the difference 
in the Lieb-Thirring inequality in two- and three-dimensions, the two-dimensional 
bound for the global attractor is a- independent. The global attractor of the two- 
dimensional LANS-a equations is similar to the bound for the global attractor of the 
two-dimensional Navier-Stokes equations. This is shown in detail in Coutand et. al. 
ETUI. 



The remainder of this chapter is organized into two sections. In Section 12.21 we 
establish the global well-posedness result, while Section I2~3l is devoted to showing the 
existence of a global H 1 attractor. 
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2.2 Global well-posedness. 

In this section, we establish the global existence of unique H l weak solutions to the 
LANS-a equations on bounded domains Q e M. 3 . Rather than using the standard 
Galerkin method, we instead take a sequence of classical solutions, via the existence 
result in [21], and prove that this sequence converges in C([0,T],iJ 1 ) to a H 1 weak 
solution of equations (j2.1.ip for all T > 0. 

2.2.1 Notation and some classical inequalities. 

We work in the following Hilbert spaces: V s = H s C\Hq and V* = {« G V s div u = 0} 
for s > 1, and V! = {m 6 V* Au = OondVl} for s > 3. We endow V* with the 
following scalar product: 

(f,9)i,a= [ f-g + a 2 Deff : Deigdx, (2.2.5) 

Jn 

where Def u is the rate of deformation tensor defined by equation 1 11. 2. 4ft . • denotes 

the usual dot product, and : is the contraction of two indices, e.g. a : b = a^bij. 

Furthermore for any integer s > 0, we set 



D s u = {D p u : \(3\ = s}, \\D s u\\ LP = ^ \\D ( 'u 



\P\=s 

where (3 — (Pi, ... , P n ) denotes a multi- index, and 



\p\ = p 1 + ... + p n , D0 = a£-~d£, 

where 

f> d 

Ok :-- 



dx k 

Throughout the dissertation, we let C > denote a generic constant. For simplic- 
ity in notation, we write u(i) = u(t, ■). The following standard inequalities are used 
frequently and we include their definitions for completeness. 
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Gagliardo-Nirenberg inequalities [TJ . Suppose 

1 i ( \ m\ , , 1 

- = - + a +(l-a)- 

p n \r n J q 

where 1/m < a < 1 (if m — n — n/r is an integer > 1, only a < 1 is allowed). Then 
for / : tt -»■ R n , 

|p7IU P <C|p m /||^-||/||ir. (2.2.6) 

Two specific cases of (12.2. 6p in dimension three are 

||^|| L 4<4||^||^ 4 ||t;||^ 4 , (2.2.7) 

WD'vWv < CWvWl-^WD^W^. (2.2.8) 

Tl 

Sobolev Embedding Theorem. If u e H s (W l ) for s > - + k, then u E C k (R n ). 

Tl 

In particular for s > -, H s (R n ) C L°°(M n ) with the bound 

||w||l°°(»") < CllwII/f'O*") (2.2.9) 

where the constant C only depends on s and n. 

Young's Inequality. Let 1 < p, q < oo, — I — = 1. Then 

p q 

ab < ea p + C(e)b p (a,b >0,e>0) (2.2.10) 

for C(e) = (ep)- q l p q-\ 

Gronwall's Inequality. Let rj(-) be a nonnegative, absolutely continuous function 

on [0, T] which satisfies for a.e. t the differential inequality 

V'(t)<<f>(t)ri(t) + il>(t), 

where <p(t) and ip(t) are nonnegative, summable functions on [0,T]. Then 



r](t) < e ^^ ds 
for all < t < T. 



rj(Q)+ I i[>(s)ds 
o 



(2.2.11) 



CHAPTER 2. Weak solutions of the LANS-a equations on bounded domains 19 

2.2.2 Three equivalent forms of the LANS-a equations. 

Three equivalent forms of the LANS-a equations will be useful to us. 
LANS-1: 

d t u + V u u + U a {u) = -vAu - (1 - a 2 A)" 1 gradp + (1 - a 2 A)" 1 / 

div u(t, x) = 0, 

(2.2.12) 

« = 0on <9f2, 
u(0, x) =u (x), 

where / G L 2 . (It is convenient to replace T in equation (12.1. lap with (1 — a 2 A)~ 1 /; 
there is no loss in generality as (1 — a 2 A) with domain H 2 n Hq is an isomorphism.) 
The Stokes operator Au = —PAu, is the Leray projection of —Am onto divergence- 
free vector fields and has domain D(A) = H 2 D Hq. As we noted above, when u = 
on dfl, then Au must also equal zero on dQ. 

LANS-2: This form is equivalent to LANS-1 in view of our remark that LANS-1 
implies Au = on dQ: 

d t (l-a 2 A)u + V u [(1 - a 2 A)u]-a 2 \7u T ■ Au = v(l-a 2 A)Au-gmdp + f, (2.2.13) 

together with the constraint div u(t,x) = and boundary data u = Au = on dfl. 
Note that when the domain Q is the period box T 3 , the Stokes operator is given by 
—A, and formulation (12.2.131) reduces to equations (12.1.41) . the LANS-a equations 
used by Foias, Holm, and Titi (T4~| . 

LANS-3: This form is the analog of the Navier-Stokes equations written in terms of 
the Helmholtz- Hodge projection: 

d t u + vAu + V a [V u u + U a {u) - (1 - a 2 A)~V] = 0, (2.2.14) 

where for s > 1, V a is the Stokes projector defined below. 
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Definition 2.2.1 For s > 1, we let V a : V s —> V* denote the Stokes projector, a 
continuous (-,-)i :a -orthogonal idempotent operator (see Proposition 1 of \%8$). It is 
defined as 

V a (w) = w — (1 — a 2 A)" 1 gradp = v, 

where (v,p) solve the Stokes problem: given w G V s , there is a unique vector field 
d6V* and a function p (unique up to an additive constant) such that 

(1 — a 2 A)v + gradp = (1 — a 2 A)w, 

div v = 0, 

v = on dQ. 

2.2.3 Results. 

We begin with two elementary lemmas. 

Lemma 2.2.1 Vf, is dense in V},. 

Proof: Let v G Vh We find a sequence v m G V 4 which converges to v in H 1 . The 
proof makes use of the Lax-Milgram Theorem to provide a compact operator from 
V^ to V 4 Define the bilinear form E : V 4 x V 4 — > R by 

E[u,v] = (u,v) H 4. 

The function space V 4 is a closed subspace of H A and therefore a Hilbert space 
endowed with the usual H A topology. Since E is defined as the usual inner product 
on H 4 , it is coercive and continuous. We define, using the Riesz- Representation 
Theorem, the bounded linear function / G ( Vt ) by 

where / G V/J. By the Lax-Milgram Theorem, there exists a unique u(f) G V 4 such 
that 

S[u,u] = (u,x;)H4 = (/,T;) ia , V?; G Vj. 
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This is the weak formulation of an elliptic problem Lu = f where L is an eighth 
order elliptic differential operator. The map L~ l : / w m is a continuous map from 
V^ — > Vt. Since H A is compactly embedded in H 1 , L~ l is a compact operator from 
V* — > Vz, and consequently has a discrete spectrum of eigenfunctions {ej}j G N- The 

oo 

eigenfunctions {ejjjgN form a Hilbert basis of V* and therefore we write v = Nc^e;. 

i=\ 

The ellipticity of L~ l implies that {ej}j g N are also eigenfunctions of L, allowing us to 

m 

conclude that {e.;} ie N G Vt- Consequently, the sequence v m = 2_] c i e i G V^ converges 

in H 1 to v as m — > oo. ■ 

||Vv||?a 
LEMMA 2.2.2 Lei Ai = inf 9 — 6e £/ie smallest eigenvalue of the Stokes opera- 

' ueV M \\ v \\h 

tor. Then for all v G V 2 , 

||Vw||i 2 + a 2 pt;||2 2 > A 1 {||t;||2 2 + a 2 ||Vw||i 2 } 

Proof: We just need to show that ||Av||| 2 > Ai||Vw|| 2 2 (since ||Vw|| 2 2 > Ai||u||^ 2 
by definition). Let {/ij}j e N be eigenfunctions associated to the Stokes operator A 

oo 

with corresponding eigenvalues Aj. For all v G V 2 , we write v = 2_. c ihi- Then Av = 



t=l 

oo 



y^CiXihi. The orthogonality property of hi implies (Av,v) = y^ c "i^i = \\^ v \\l 2 an d 

i=l i=l 

consequently, 

oo / oo \ 

i=l \t=l J 

which completes the proof. | 

Weak solutions to the isotropic LANS-a equations ( 12.1.11) are defined below. 

Definition 2.2.2 Let f e L 2 (fi) and u G V*. For any T > 0, a function 

ueC([0,T};Vl)nL\[0,T];D(A)) 

with 4s G L 2 ([0,T];L 2 ) n L°°((0,T); V*) and w(0) = u zs sazd to 6e a wea£ solution 
to the LANS-a equations with initial data uq in the interval [0,T] provided 
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-(l-a 2 A)u,w) + u((l-a 2 A)Au,w) 

at I D{A) 

+ (B(u,(l-a 2 A)u),w) D{A) = (f,w) (2.2.15) 

for every w = w(x) G D(A), the domain of the Stokes operator, and for almost every 
t G [0, T] with 

B(u,v) = V u v + Vu T -v] 

moreover, u(0) = u in V*. Here, the equation Ii2.2.15\) is understood in the following 
sense: for every t ,t G [0,T], 

(u(t), (1 - a 2 A)w) D(A) + J (B(u(s), (1 - a 2 A)u(s)),w) D{A) ds 

= - u(Au(s),(l -a 2 A)w) D(A) ds+ / (f,w)ds. 

J to J to 

For our proof, we shall make use of the following results from Marsden and Shkoller 

EH 



Lemma 2.2.3 (Theorem 5.2 of |2I]) For u e V", s G [3,5), and f e L 2 (Q), 
there exists a unique solution u to equation Ii2.2.12\) in C([0, oo), V*). 

and 

Lemma 2.2.4 (Lemma 5.1 of |2T]) For s > 3, U a : V s -> V s andU a : V 2 -»■ F 1 ^ 
forae (0,|). 

We can now state our main result. 

Theorem 2.2.1 For f G L 2 (f2) and Uq G V^, i/iere exists a unique weak solution 

u G C([0,oo), Vi) nL°°(( S ,oo), Vj), Vs > 
to equation $2.2. 12\) . The solution depends continuously on the initial data uq. 
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Proof: Consider a sequence of initial velocity {u e Q }% l G V^ and force {f € }% i G i? 1 
such that Uq — > uq (by Lemma [2.2.11) and f e — > /. By Lemma [2.2.31 for each e G N 
there exists a unique solution w e of the LANS-a equations in C([0, oo), V*). This 
solution has sufficient regularity to satisfy the weak formulation: 

<?/(£), (1 - a 2 A)w) + [ (B(u*(s), (1 - a 2 A)u £ (s)), w) ds 



to 



- / v(Au e (s),(l-a 2 A)w) ds + / (f e ,w)ds 



D(A) 
to J to 

for all w G D(A). In order to show the limit as e — ► oo is also a weak solution, we 
need to develop the appropriate energy estimates. 

An H l Estimate. Define a € = V u ^u e + U a (u € ) - (1 - a 2 A)- l f\ Since u e G H% we 
may use it as a test function. Taking the (•, -)i iQ inner product, defined by ( 12.2.51) . of 
LANS-3 together with u e , 
d 



—u e ,u £ 



dt 

By the definition of (•, -)i iQ , 
d 



+ u{Au\u e )^ a + (P Q (a e )X)i,a = 0. 



l,a 



1 d 



u u 6 ,u e ) +z/(An e ,^)i, a ==^[||n e ||i 2 + a 2 ||Vn e ||i 2 ]+z/[||V^||i 2 +a 2 ||An e ||i 2 ] 

Notice that the previous integration by parts used the additional boundary condition 
Au e = 0. Using the properties of V a , we find 

(P a (a%v*) lia = -(f e ,u e )- 



Therefore, 

\\u e \\ 2 L 2 + a 2 \\Vu e \\ 2 L2 ] +z/[||VM e || 2 :2 +a 2 ||AM e || 2 2 (f,u e 
By Poincare's inequality and Young's inequality ( 12.2. 101) . 



ld_ 
2dt 



(2.2.16) 



d_ 

dt 



.e I |2 



u6 \\h + a l|VM e ||i2 +u ||Vm £ ||22 + a z \\Au £ \\ 2 L 



.e| (2 



.el 12 



<K t 



(2.2.17) 
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where K\ 



lint 



v\ 



&-. By Lemma EJL21 



d 
db 



I ei|2 | 2i|v7 ei|2 
\U || L 2 + a || Vlt 1 1 L 2 



+ ^Ai 



I ei|2 , 2i|v7 e 1 1 2 
P IIl 2 + a II VU ||x,2 



Using Gronwall's inequality (j2.2.1ip . we get the if 1 estimate 



<#i. 



#i 



\u%t)\\ 2 L2 + a 2 \\Vu%t)\\ 2 L2 < e^ i (||^(0)||i 2 + a 2 ||V^(0)||i2)+-^(l 

< h:=\K\\ 2 L , + a 2 \\vu^o)\\^ + 



-i/Ai t> 



uXi 



(2.2.18) 



where Ai is defined in Lemma [2.2.21 Therefore we have proved that 

u e eL~([0,oo),Vj) 

independently of e. 

An H 2 Estimate. Since Au e is divergence-free, we take the L 2 inner product of 

LANS-2 together with Au e to get 

1 d 



Idt 



|Vw e ||| 2 + a 2 \\Au e \\ 2 L2 +u[\\Au e \\ 2 L 2 + a 2 \\VAu e \\ 2 L 2\ 

+ (B(u e , (1 - a 2 A)u £ ), Au e ) = (f e , Au e ). (2.2.19) 



By estimate (12.2. 18j) . we use the Sobolev embedding result W l, °° C L 4 to bound the 
nonlinear term of equation (12.2.191) 

|(V u ,u e ,Au e )| < ||W|| L 2||u e || L 4||y4u e || L 4 < CA;i||ii 6 || H i||Au e || L 4 



1 3/4 1 



1/1 



< Ck{\\VAu e \\22\\Au e \\i; 

< ^!|| V ^||i 2 + ^||^||| 2 + c(iKiUx), 



(2.2.20) 



may be estimated by the Gagliardo-Nirenberg inequality (12.2.61) and a repeated use of 
Young's inequality. The boundary conditions, the incompressibility constraint, and 
one integration by parts lead to 



\{V u ^u\Atf)\=U. 



(2.2.21) 



CHAPTER 2. Weak solutions of the LANS-a equations on bounded domains 



25 



For the third term, using (12.2.61) . (12.2.71) . and (I2.2.18P we have the estimate, 



((W T )- Au e ,Au e ) | < \\VAu e \\ L 2\\u € \\ L 4\\Au e \\ L 4<C\\u e \\ H i\\VAu e 



i7/4, 

\L 2 I 



\Au e 



,1/4 

II-' 



< CkiWVAu 6 



1 23/12 1 1 eiil/12 



lL-> 



II 



,13/12i 



\j2 _ ^ *V1 V ^T-li I I r 2 



23/12 



V, 



-l|v^||| 2 + c(||« e || H1 ; 



(2.2.22) 



Again Young's inequality implies 



\(r,Au<)\ < nriu 2 |i^n L2 < l\\r\\h + ^ii^iii- 



Therefore by (I2.2.20I) - (I2.2.22I) . the equation 02.2.19P becomes 



1 d 

2di 



Vu e \\ z L 2 + a 2 \\Au 



e 1 1 2 

II 2 



is [\\Au e \\ 2 L 2 + a 2 \\VAu £ \\ 2 L2 ] 



< \(B(u',(l-a 2 A)u*),Au*)\ + \(r,Au')\ 

< ^\\vAu<\\i, + l\\Avr\\b + h\nh + c. 



Hence by Lemma 12.2.21 we have 



d 
di 



\Vu e \\ 2 L 2 + a 2 \\Au 
d 



£||2 
L 2 



+ v\i 



V 



e I 1 2 



W||£ 2 +a z \\Au e \\l 2 



<K? 



(2.2.23) 



- [(||Vn(t)||i 2 + a 2 ||^(t)||i 2 )e^] <K 2 e^. 



where K 2 = ^||/ e ||^ 2 + C. Smce\\Au\\ 2 L2 is not necessarily bounded at t = 0, for every 
Si > we integrate both sides from [si,£], 

liv^(t)ni 2 + « 2 n^^(t)||| 2 < (l|v^( Sl )lli 2 + « 2 n^^( Sl )lli 2 )e" Al(si " i) + ^a x ^ 2 . 

(2.2.24) 
Thus, independently of e, 



^GL°°(( Sl ,oo),V 2 (fi)), 

for all s\ > 

An H 3 Estimate. Since A 2 u e is not necessarily equal to zero on dQ, we do not use it 

to derive an i/ 3 -estimate. To achieve this estimate we make use of the Ladyzhenskaya 
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method [18] . Let u t denote d t u, and differentiate the equations LANS-2 with respect 
to time to get 

d t (l-a 2 A)u e t + V ul {l-a 2 A)u e + V u .{l-a 2 A)u e t -a 2 {Vu e t ) T ■ Au e 
- a 2 (Vu e t ) T ■ Au\ = -gradp t - u{\ - a 2 A)Au e t . 

Noting that u\ E D(A), we take the L 2 inner product with u t to get 

\j t l\\<Wl> + oc 2 \\Vut\\l 2 ]+u[\\Vut\\ 2 L2 + a 2 \\Aui\\l 2 ] 

< |(^V^)|+a 2 |(A<,V^ + V M ^)|. (2.2.25) 

Before estimating the right hand side, we would like a bound for 1 1 u|| j ^ 2 which will 
be used in the later computations. For each e, u e satisfies the equation LANS-3 and 
consequently, 

|K e ||| 2 < u\(Au^ul)\ + \(V u ,u%ul)\ 

+ |(W> e )X)| + |((1 - a 2 AY l r ,ut)\. (2.2.26) 

The first term is simply v^Atf , uf)\ < |||w|||^ a + 2^ 2 ||Au e || 2 : 2 by Young's inequality. 
Using (I2.2.7P and Young's inequality, the second term becomes, 

|(V u eU e X)| < 16||Vn £ || L2 ||n e ||^ 4 ||^|| L2 pn e ||g 4 

< v 2 \\^\\ 2 L i + l\\ul\\ 2 L2 + Ck{ , 

where k\ is the time independent H 1 bound (12.2.181) . Let v e = (1 — a 2 A)~ l u\. Then, 
by integration by parts, the third term of (I2.2.26P is 

\{U a {u e ),u\)\ < a 2 \(Au e -Vu e ,v e )\<a 2 \\u € \\ 2 H2 \\v € \\ H i 



a 



2 



,e 1 12 i „,2||„.e||4 



, u- ilff 1 T" II" Hi? 2 " 
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By definition, v e — a 2 Av e = u\. Taking the L 2 inner product of both sides with v t we 
have an H 1 estimate 

M\h + « 2 liw||i 2 = «x) < \\\ut\\l> + \\\v%l 2 . 

Hence ||w e ||^i < ^ll^illi 2 and therefore 

\{U a ^)X t )\<\\\ul\\h + a 2 \\u^\\ A H2 . 
The last term of (12.2.261) can be estimated in a similar way 

|((l-a 2 A)- 1 A<)l = l(A^)l<^IKIIi 2 +2||r||i 2 <^||^||i 2 + 2||r||i 2 . 
Combining the above estimates, we obtain the bound 

| \u e t \\ 2 L2 <Qv 2 \\Au e \\l 2 +2a 2 \ |u £ |& 2 +4||f ||| 2 + Cfci°, (2.2.27) 

From the H 2 estimate (12.2.24J) we conclude that u\ G L°° ((si, oo), L 2 ) for all s 1 > 0. 
Using this result, we now estimate each of the terms on the right hand side of (12.2.25)) . 
Using (I2.2.6P - 02.2.8I) and Young's inequality, we get 

|«,V u| OI + a 2 \(AulV u eU e t + V ut u e )\ 



< C 



V7 e 1 1 3/4 , , e||5/4|| ell 2|| A el 1 15/8 1 1 e||l/8|i E || 

V«t|| L a |pt|| L a \\U \\ H 2 + a \\Au t \\ L i \\u\\\£ 3 \\U \\ H i 

3/4 1 , ,ei|l/4| 



i 211 4 ell llvr ellJ/^ll enV^II ell 

+ a \\Au t \\ L 2\\\7u t \\£ 2 ||u || L 2 ||w ||ffi 

2 

< ^11^7,^1 1 2 4- Va II47/ £ II 2 4- r 
^ 2" *" L _ 9 _ " *" L Sl ' 

where C Sl is a constant multiple of the H 2 bound depending only on the lower time 
bound S\ > 0. Therefore (12.2.251) becomes 

By Lemma E23 

^[INlL> + ^||Vt4||y+i/A 1 [||t4||i a + a a ||Vu|||ia]<C7. 1 
^[(|K||l 2 + a 2 ||V^||i 2 )e^]<C Sl e^. 
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For every s 2 > 0, we integrate from [s2,t] and obtain the H 1 bound of u|, 

ll^(t)ll| 2 + « 2 ||v^(t)||| 2 < (ii^c^^nl. + ^iiv^c^)!!^)^ 1 ^-*^^!^, 

Hence, independently of e, 

u t e eL 00 ((s 2 ,oo),Vj), 

for all s 2 > 0. 

Let s = min{si, ^j- By standard Sobolev inequalities, we also have V u ^u e belongs 
L°°((s,oo),iJ 1 ). Therefore, 

\\V a {V u .u e )\\ v i < \\V u .u € \\ m , 

implies that for all s > 0, V a (V u ^u t ) belongs to L°°((s, 00), V*), independently of 
e. From Lemma 12.2.41 and our H 2 estimate (j2.2.24p . we infer that U a (u e ) is con- 
tained in L°°((s, 00), H 1+cr ) for < a < 1/3 and consequently V a (U a (u £ )) belongs 
L°°((s, 00), V* +cr ), independently of e. 

Thus using the equation LANS-3 and the results for u|, we have that for any 
s > 0, 

i/A« e = -V a [V u ,u e +U a {u e ) - (1 - a 2 A)" 1 / 6 ] + u\ 

belongs to L°°((s, 00), V*), independently of e. By the elliptic regularity of the Stokes 
operator A, we conclude that 

u e eL~((a,oo),iS). 

Define # = {w e L 2 | div w = 0}. For all T > 0, 02.2.271) implies u\ e L 2 ([0, I 1 ], #) 
(independently of e) and therefore the classical compactness theorem (see for in- 
stance [131 HH]) enables us to conclude that there is a subsequence u e such that for 
all T > 0, 

u £ ' -»■ u weakly in L 2 ([0,T],D(A)), 
u 6 ' -► u strongly in L 2 ([0, T], V*), 
m £ ' -»• u in C([0,T],if). 
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It is straightforward to verify that u satisfies the weak formulation associated with 
/. Indeed, from an identical argument provided in [13], we may conclude that u G 
C([0, oo), V*). Furthermore, we infer that u(0) = Uq. Hence we are lead to the 
existence of weak solutions for the LANS-a equation. The uniqueness and continuous 
dependence on initial data of weak solutions can be proved in the same classical way 
as was done in [13] on Page 14 for the periodic case. | 

Remark 2.2.1 A standard contraction mapping argument can be used to show the 
existence of LANS-a solutions u in C([0, T],V*) for s > 1. In the absence of forcing 
or when f is C°° , the solution u of the LANS-a equations is instantly regularized so 
thatu e C°°((0,oo) x ft). 

2.3 Estimating the dimension of the global attractor. 

Theorem 12.2.11 is sufficient to define the semi-group S{t) by 

S{t) : n e v; h+ u (t) e v;. 

We now show the uniform compactness property of the operator S and the existence 
of an H 2 absorbing set. By Theorem 1. 1.1 in [29], this implies the existence of a 
maximal compact global H 1 attractor. 

2.3.1 Absorbing sets and attractor. 

In this section we prove the existence of an absorbing set in V 2 . Usually proving 
the existence of absorbing sets amounts to proving a priori estimates; for the LANS- 
a equations, these estimates are established in Section 2.3. Thus estimate (12.2 .241) 
implies 

limsupca 2 |H|^2 < p\ := u\ 1 K 2 = 2Ai||/|||a + C(\\u\\ H i). 

t— >+oo 
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We conclude that the ball -By 2 (0, p 2 ) of V 2 , denoted B 2 , is an absorbing set in V 2 for 
the semi-group S(t). We choose p 2 > P2 an d denote B the ball B v i(0, p 2 ). If & is 
any bounded set of V*, then from (12.2.241) . S{t)B C £ for t > t (B, p 2 '), where 

t = Si + — T- log J5 2 • 

1/A1 p 2 ' - P2 

At the same time, this result proves the uniform compactness of S(t); any bounded 
set B of V* is included in such a ball, and for w £ V* and t >t , u(t) belongs to £>2 
which is bounded in H 2 and relatively compact in H 1 . Consequently, by Theorem 

1. 1.1 in [29], the LANS-a equations have a nonempty compact, convex, and connected 
global attractor in H 1 , 

A a = n s>0 (u t > s S(t)B 2 ). 

2.3.2 Dimension of the attractor. 

The dimension of the three-dimensional attractor for the LANS-a equations follows 
the arguments of Theorem 6 in [2] (stated below), and leads to the same bound 
as in the periodic case. It is important to note that the bound tends to infinity as 
a^O. 

Theorem 2.3.1 (Theorem 6 of [T3]) The Hausdorff and fractal dimensions of 
the global attractor of the LANS-a equations, dn{A a ) and dp{A a ), respectively, sat- 
isfy: 

MA') < MA") < <w{g* (^fV 2 (^ f ) 3 ' 



I f\ 

where G = ^ 4 is the Grashoff number, and d > is a constant depending only of 
the shape ofVt. 

In two dimensions, [10] proved that the dimension of the attractor can be bounded 
independently of a by the Grashoff number. Hence, in two- dimensions, we find a 
similar bound as for the Navier-Stokes equations. 
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Theorem 2.3.2 (Theorem 3 of [TU]) The Hausdorff and fractal dimensions of 
the global attractor of the LANS-a equations, dff(A) and dp(A), respectively, sat- 
isfy: 

d H (A) < d F (A) < C'G, 



where G = ^-^- — is the Grashoff number, and C > is a constant depending only 
on the shape offl. 
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Chapter 3 

Local well-posedness of the 
anisotropic LANS-a equations on 

T 3 



3.1 The anisotropic LANS-a equations. 

The anisotropic Lagrangian averaged Navier-Stokes (LANS-a) equations are given on 
Q x (0, T) by 

(1 -a 2 C) (d t u + div u <g> u - uFCu) = -gradp, (3.1.1a) 

divu(£,x) = 0, (3.1.1b) 

d t F + VF-u = Vu-F+ [Vu-Ff, (3.1.1c) 

m(0, x) = u (x), F(0, x) = F {x) > 0, (3.1. Id) 

where u(t,x) denotes the divergence-free mean velocity vector, p(t,x) represents the 
scalar pressure field, v is the kinematic viscosity, F(t, x) denotes the covariance tensor 
(a 3 x 3 matrix), and the notation Fq > means the initial covariance matrix is 
assumed positive semi-definite. In [22], the authors define the covariance tensor F as 
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the ensemble average of the tensor product of the Lagrangian fluctuation vector £' 
given explicitly as 

F(t,x) = (C(t,x)®e(t,x)). 

Furthermore, Cu := div[Vw • F] and the operator P is a projection onto divergence- 
free vector fields. When F is assumed to be a constant multiple of the identity matrix, 
the anisotropic equations reduce to the isotropic equations of Chapter [2] for only the 
averaged fluid velocity. A brief history of the isotropic equations can be found in 
Chapter El 

The operator P can be chosen to be either the Leray projector P : L 2 {Q) — > 
{v G L 2 (Q) | divx; = 0, v ■ n = on dQ}, or the generalized Stokes projector Vp := 
[P(l — a 2 C)]~ 1 P(l — a 2 C) defined in detail by the following definition. 

Definition 3.1.1 For any v in the domain of (1 — a 2 C) and F > 0, we set the 
generalized Stokes projector Vp(v) = w, where w is the solution of the generalized 
Stokes problem 

(1 - a 2 C)w + gradp = (1 - a 2 C)v 
div w = w = on dQ; 

thus 

p°( v ) = v - (1 - a 2 C)^ 1 gradp. 

Remark 3.1.1 The Stokes projector defined in Definition \2.2.1\ is a special case of 
Definition \3.1.1\ with the condition F = Id. 

The generalized Stokes projector Vp maps onto the space of H l divergence-free 
vector fields that vanish on dQ. Since the Lagrangian fluctuations are necessarily zero 
along dQ, we know F(t,x) = on dQ for all time t > 0. Unlike the Leray projector 
which only assigns the normal component of u to vanish on the boundary, the gener- 
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alized Stokes projector provides the correct boundary condition for F. Therefore on 
bounded domains, it is important to use this projector. 

Modeling incompressible flow on bounded domains is the purpose of the anisotropic 
model. The first step towards this goal is proving well-posedness of the equations on 
unbounded domains. In this chapter, we shall restrict our attention to the three- 
dimensional periodic box Q = T 3 . We define for s > 0, 

H s pcr = {u E [H S (T 3 )] 3 | f u(x)dx = 0} 

Jr 3 

and 

H s div = {ueH; cr \ divu = 0}. 

Marsden and Shkoller [22] prove the local well-posedness of the anisotropic LANS-a 
equations with the projector P = P for initial velocity fields in the class H div and 
positive F G [ifp Cr (T 3 )] 3x3 for s > 3.5. They use the classical Galerkin projection 
to obtain smooth finite-dimensional approximations u m and F m . After establishing a 
priori bounds independent of m, the authors appeal to compactness results to extract 
a subsequence of u m which converges in C°([0, T]; i?| iv ) uniformly in m. An important 
tool in obtaining the proper estimates is the control of the regularity of the covariance 
tensor F by the one derivative higher regularity of the mean velocity u. This allows 
the authors to establish estimates only dependent on the initial covariance tensor. 

The next step towards applying the anisotropic equations to bounded domains 
is the study of solutions to the equations with the generalized Stokes projector vis- 
cosity term. As discussed above, the generalized Stokes projector preserved the zero 
boundary condition for the covariance tensor. Although we known the correct projec- 
tion for bounded domains, the appropriate additional boundary condition necessary 
to solve the forth order equation A3. 1. lj) is unknown. In this chapter, we show the 
well-posedness of the anisotropic LANS-a equations on T 3 with the projector P = Vp 
by a method similar to the one used by Marsden and Shkoller [22]. However, in ob- 
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taining a priori P s -type bounds, it is important that we use the differential operator 
P°P(1 — a 2 C), \a\ < s — 2, rather than the standard D s , where D s denotes all partials 
derivatives of order s. This particular approach is different than what is done in |22j . 
It allows us to estimate the viscosity term without commuting the operator D s with 
the inverse operator [P(l — a^C)]" 1 and leads us to bounds involving the Sobolev 
norms of P(l — a 2 C)u and u. We then appeal to the elliptic regularity estimates 
established in Proposition 13.2. II to write all bounds in terms of P(l — a 2 C)u. The 
same compactness argument leads us to the existence of solutions to the anisotropic 
LANS-a equations ( 13.1.11) . 

3.2 Local well-posedness with Stokes projector. 

Before proving the existence and uniqueness of classical solutions to the anisotropic 
LANS-a equations stated above, we establish an elliptic regularity estimate for the 
operator P(l — a 2 C) which holds assuming the least regularity of P. The proof relies 
on a standard commutator estimate which we include below for completeness. For 
the remainder of the chapter, C represents an arbitrary positive constant unless, by 
the context, it is clear that C is the operator defined above. 

Lemma 3.2.1 For f,g e H s (n), 

WlD s ,f]g\\ L 2<C(\\f\\H°\\g\\L™ + ||V/|U=c|b||H.-i) (3.2.2) 

Proposition 3.2.1 For P > and F e [Pp" 1 ^ 3 )] K , suppose u e H k (T 3 ) solves 
P(l - a 2 C)u = w fork>3. If w e H^T 3 ), then 

||m|Ih3<C , (1 + C , ||P||]| 2 )||w||^ 1 , (3.2.3) 

andforw eH k ~ 2 (T 3 ), 

WuWl^CWwWl^ + CWFWl^WuWl,.^ k>4. (3.2.4) 
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Moreover, if F solves 113.1.1]) then for some T > 0, 

||*3 < C^HFoll^lHllri + Cdl^oll^lHlHi, 
\\u\\ 2 m < C(tA\F \\ H 2)\\w\\ 2 H , + C(\\F \\ H3 )\\w\\ 4 H1 , 



(3.2.5) 
(3.2.6) 



\u\\%k < Gllwll^-a + C^llFoll^-ilHI^-x + CHFoll^-ilHI^-i, (3.2.7) 



for k > 4 and all t <T. 

Remark 3.2.1 In a classical elliptic regularity statement one assumes u G H 1 solves 
the elliptic second order differential equation and attempts to prove u has two deriva- 
tives more regularity than the inhomogeneous term using difference quotients. Since 
we are assuming that ||F||x,oo is not necessarily bounded, this regularity is not guar- 
anteed and we add the assumption u G H k . Later in the chapter, the proposition is 
used when u G C°°([0, T] x T 3 ) and therefore the regularity of u is not delicate. 

Remark 3.2.2 If F e [C°°(T 3 )] 3x3 , then equations (EJLlj) and (3~J^\) reduce to the 
standard elliptic regularity result 



\\u\\ 2 H k < C(||F|| W fc-i,oc»)||w||^fc_i. 
Proof. We first write P(l — a 2 C)u = w in the energy form 



(3.2.8) 



{u, v) L2 + a 2 {F ■ Vu, Wv) L 2 = (w, v) L 2, v G # p fe er (T 3 ) k > 1. 



(3.2.9) 



We begin by showing the special case of k — 3, namely inequality (13.2.31) . When 
v = u, uniform ellipticity of F implies (F • Vu, Vu)z,2 > A||u||^i from which we find 



\\u\\ 2 H i <C|H|#_i. (3.2.10) 

Now suppose v = D 4 u and input this into inequality (J3.2.9P - Integrating by parts 



gives 



a\D\F ■ Vu), D 2 Vu) L 2 < (w, D%) L 2. 
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Expanding the left hand side, 

a 2 (D 2 (F-Vu),D 2 Vu) L 2 = a 2 (F ■ D 2 Vu, D 2 Vu) L 2 + 2a 2 (VF ■ DVu, D 2 Vu) L 2 

+a 2 (D 2 F ■ Vu, D 2 Vu) L 2 = : / + II + ///, 

where / > a 2 A||M||^ 3 by uniform ellipticity of F and 

II < C||VF|| L 4||L) 2 n|| L 4||D 3 u|| L 2 < C||F||^2||D 2 n|| L 4||n||^, 

III < C7||Z> 2 i^||i3||Vu|jioo||Z) 3 M|| £ a <C||F|| H 2||u|| W 2,4||u||if3 
< C||F| \ H 2 | |D 2 ti|| L 4 ||w||#3, 

by standard Sobolev inequalities. Combining the estimates above, we get 



H\ 2 m 



< C(w, D*u) + C\\F\\ H 2\\D 2 u\\ L 4\\u\\ H 3 

^ /~i\\ II II 7^»4 ii i /nri I tttI I || || 7/8 1 1 1 1 1/8 1 1 1 1 

< G||w||#i||.L> u\\ H -i +G||i < ||_ ff 2||ti||^3||ii||^ 1 ||u||_ H -3 

< £\H\ 2 HS + c\\w\\ 2 H1 + cwfw^WuW 2 ^, 

where we have used the Gagliardo-Nirenberg inequality (I2.2.6P for the second inequal- 
ity and Young's inequality (12.2.101) for the last. Using estimate (13.2.101) and choosing 
e > sufficiently small, we can absorb the H-uH^-term into the left side, obtaining 

||«|||.<t7(l + C7||F||^)|HI|i. 

For k > 4, we substitute v = (—l) k ~ l D 2k ~ 2 u into (I3.2.9p . After integrating by 
parts, 

a 2 (D k ~\F-Wu),D k - l Wu) L 2 < (D k ~ l w,D k ~ l u) L 2. 

Expanding the left hand side we see 

a 2 {D k ~ 1 (F-Vu),D k ~ 1 Vu) L 2 = a 2 {F ■ D k ~ 1 Vu,D k - 1 Vu) L 2 

+a 2 ([D k -\ F]Vu, D k ~ l Vu) L 2 =: I + II, 
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where / > a 2 A||n||^. fe . Unlike the case when k = 3, /c > 4 > | implies H k ~ l C W 1 ' 00 
and therefore by a standard commutator estimate (13.2.21) . 

II < C (\\F\\ H h-i\\\I U\\l°° + ||Vi r '||z / °o||VM||fl-fc-2)||u||fl'fc < C||i 7 '||_ H -fc-i||w|| fl 'fc-i||u||_ ff jfe. 

Hence 

\\u\\ 2 H k < C\\D k ~ 1 w\\ H -i\\D k ~ 1 u\\ H i + C\\F\\ H k-i\\u\\ H k-i\\u\\ H h, 
< e\\u\\ 2 Hk + C\\w\\ 2 Hk -2 + C\\F\\ 2 H k-i\\u\\ 2 Hk - 



fe-i. 



By choosing e > sufficiently small we achieve inequality (13.2.41) . 

To show inequality (I3.2.5p . recall from [22], equation (33), that if F solves (13.1. Id) 
for a given u, we may estimate ||.F||#s by one higher order derivative of u in the 
following way, 

\\F\\h* < C||Fo|| H .exp / \\u(s)\\ H s+i ds < C\\F \\ H s (l + t\\u\\ H s+i + 0(t 2 )) . 

Jo 

(3.2.11) 

For k — 3, it follows that inequality ( I3.2.3P can be written as 



Mltf3 



2 m < c\\w\\ 2 H1 + c\\F \\^ 2 (i + m\\u\\ H3 + o(t 2 ))\\w\\ 2 Hv 

< C\\w\\ 2 Hl +C\\F \\^ 2 (l + O(t 2 ))\\w\\ 2 Hl + tC\\F \\ 1 £ a \\w\\ 2 H x\\u\\ H 3, 

< C\\w\\ 2 Hl + C(t)\\F \\™\\w\\ 2 Hl +C\\F \\ 32 2 \\w\\% 1 + et 2 \\u\\ 2 m , 

For t < T, we choose e > sufficiently small to conclude inequality ( 13.2.5J) . The 
estimate for k > 4 follows from the same procedure. I 

Theorem 3.2.1 For s > 7/2, and u e H^ iv (T 3 ), F e [^ cr (T 3 )] 3x3 , with F > 0, 
there exists a unique solution (u, F) with u £ C°([0,T]; ifj iv ) fl L 2 (0,T; H^ 1 ) and 
F G C°([0,T]; [-ffp Cr ] 3x3 ) to the anisotropic LANS-a equations, where T depends on 
the initial data. 
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Proof. The evolution equation for the mean velocity (13.1. lap may be written as 

d t u + V^ div(w <g> u) = vV^Cu (3.2.12) 

u(0,x) = u (x). 

Approximate solutions. Let v,i = V{(-, x) (i = 1,2, . . .) denote the smooth periodic 
orthogonal basis of H^ CT given by the eigenfunctions of the Stokes operator. Define 
P m : #p Cr -> V m := span{ui, . . . , v m } and u m = P m u. 
Consider the Galerkin projection of (j3.1.ip given by 

d t u m + PmPp m div(w m <g> u m ) = vP m Vp m C m u m , 

(3.2.13) 

u m (0,x) = P m u , 
with C m u m := div(F m • Vu m ) where F m is a solution to 

(3.2.14) 
F m (0,x)=F (t,x)>0. 

For each meN, there is a smooth solution; in order to pass to the limit as m — > oo 
to produce a solution to (13.1.11) . it suffices to obtain a priori estimates independent 
of m. 

H s estimate. Since each term on the right-hand side of equation (13.2. 13p involves 
the generalized Stokes projector Vp m := [P(l — a 2 C m )] P(l — a 2 C m ), we use the 
operator P(l — a 2 C m ) before the standard differential operator in the following H s - 
energy estimate. By letting a denote a multi-index with \a\ < s — 2, we apply 
D a P{l — a 2 C rn ) to both sides of equation (13.2.131) and take the L 2 -inner product with 
D a P{l — a 2 C m )u m . The left-hand side reduces to 

{D°P{1 - a 2 C m )d t u m , D°w m ) L 2 = (D°d t w m , D°w m ) L 2 - a 2 (D°P[C m , d t ]u m , D°w m ) L i 

= {d t D°w m ,D°w m ) L 2 + a 2 (D (7 Pdiv(d t F m ■ Vu m ),D°w m ) L2 
= ~ t \\ Dffw rn\\h-a 2 (D< T P(d t F m -Vu m ),VD ff w m ) L 2, 
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where we define w m := P(l — a 2 C m )u m . In addition 

j t \\D°w m \\ 2 L2 = 2a 2 (D°P(d t F m -Vu m ),VD°w m ) L 2 

-2{D°P{1 - a 2 C m ) div(« m ® Um), D°W m ) L 2 
+2v{D°P{l - a 2 C m )C m u m ,D°w m ) L 2 
= I + 11 + III, (3.2.15) 

where we have used the definition of Vp m := [P(l — a 2 C m )] _1 P(l — a 2 C m ) for the II 
and /// terms. 

Estimate for I. By the continuity of the Leray projector P and the fact that H r is 
an multiplicative algebra for r > |, we use equation (13.2. 14p to write the first term 

as 

\I\ < C(a 2 )\\P(d t F m - Vu m )\\ H .-2\\w m \\ H .-i < C(a 2 )\\d t F m \\ H s-2\\u m \\ H s-i\\w m \\ HS -i 

< C(a 2 )(\\VF m ■ u m \\ H s-2 + ||Vw m • i r m || ir »-2)||« m || H - 8 -i||iy m | 1^,-1 

< C(a 2 )\\F m \\ H s-i\\u m \\ 2 H s-i\\w m \\ H s-^ 

< £i\\w m \\ 2 H s-i + C(e 1 ,a i )\\F m \\ 2 H s-.i\\u m \\ i H s-i, 

where we have used Young's inequality for the last estimate. 
Estimate for II. Since div D°~w m = 0, 

\II\ < C\(D s - 2 (l-a 2 C m )dw(u m ®u m ),D s - 2 w m )\ 

< C\\<&v(u m ®u rn ,)\\H*-2\\w m \\ H *-* 
+C(a 2 )\(D s - 2 (F m -Vdiv(u m ®u m )),VD s - 2 w m )\ 

< C\\Vu m ■ u m \\ Hs -2\\w m \\ H s-2 +C(a 2 )\\F m ■ Vdiv(w m ® u m ))\\ H s-2\\w m \\ H s-i. 

By Young's inequality we conclude 

\II\ < e2\\w rn \\ H s-i+C\\u m \\ H s-2\\u m \\ H s-i\\w m \\ H s-2+C(a )\\F m \\ Ha -3\\u m \\ H ,-i\\u m \\ Ha 

Estimate for III. Using the coercivity of F m , the estimate for expression III includes 
a negative coefficient in front of the ||iu m ||#s-i-term. The negative coefficient allows 
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us to control the ||io m ||# s -i-terms from estimates for expressions I and II. 

/// = 2v{D°P{l - a 2 C m )C m u m , D°w m ) L 2 = 2v{D°PC m {l - a 2 C m )u m , D°w m ) L . 
= 2v(D°C m P(l - a 2 C m )u mi D°w m ) L 2 + 2v(D°[P, C m ](l - a 2 C m )u m , D a w m ) L 2 
-2v{D°{F m ■ Vw m ), VD a w m ) L 2 + 2v(D°[P, C m ]u m , D a w m ) L 2 
-2a 2 v{D°\P, C m ]C m u m , D a w) L 2. (3.2.16) 

Since d\vu m = 0, the definition of the Leray projector P implies [P, C m ]u m = 
PC m u m — C m Pu m = — gradpi where Apx = div C m u m . Therefore the second term 
on the right hand side of equation (13.2. 16j) vanishes. Also 

[Pj ^m\^mU m = P^m U m ~ ^mP^m u m = ^ m u m ~~ grad£>2 ~~ ^rn U m + ^m . grad^i 

= - gradp 2 + gradC m j»i + [C m , V]pi =: gradP - div(VF m • Vpi), 
where Ap2 = divC^w m . Using Leibniz' formula on the first term of (13.2.161) we find 
/// = -2v{F m ■ VD°w m , VD°w m ) L 2 + C ^ (D p F m D^Vw m ,VD a w m ) L 2 

0<|/3|<k| 

-2a 2 v{D°{VF m ■ Vpi), VD a w m ) L 2. (3.2.17) 

Rather than introducing the commutator in equation ( 13. 2. 17ft . we use the Leibniz' 
formula to obtain the correct power on the term ||w m ||_H"3 for the critical case of 
a = 4. With a standard commutator estimate the best we can obtain is a bound by 
||-Fm||.ff 3 ||wm||fl-3 ) which doesn't allow the correct isolation of the ||iu m ||#3-terin. 
By the uniform ellipticity of F m , the first term of ( 13.2. 17ft 

-2u{F m ■ VD a w m , VD a w m ) L 2 < -\\\w m \\ 2 H .-i. 
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The second term of inequality (J3.2.17P 
Y^ {DPF m D°-PVw m ,VD° Wm ) L *<C J2 \\D p F m \\ L *\\D°-f> +1 w m \\ L *\\VD''w m \\v 

0<|/3|<H 0<|/3|<|<r| 

< C\ \D s ~ 2 F m \ \ H i | \D s ~ 2 w m \ \ L 4 1 \w m \ \ H .-i 

3 1 

< C\\F m \\ H s-i\\D s ~ 1 w m \\l 2 \\D s ~ 2 w m \\l 2 \\w m \\ H s-i 



< C\\F m \\ Hs -i\\w m \\^.,_ 



7 
s—2 1 1 UJ m\ \jjs—l i 



where we have used the Gagliardo-Nirenberg inequality for the third inequality. Fi- 
nally the last term of (J3.2.17P is 

-2a 2 u(D a (VF m -Vp 1 ),VD' T w m ) L 2 < C(a 2 )\\F m \\ H s-i\\V P i\\ H s-2\\w m \\ H s-i 

< C{a )\\F m \\ H s~i\\C rn u m \\ H s~2\\w rn \\ H s-i 

< C(a )||-F m || Hs _i||'u m ||H»|l w m||H s - 1 

Hence by Young's inequality we combine the above estimates to write inequality 
03.2. 17j> as 

\III\ < (-A + e 3 )||w m ||^_ 1 +C||F m ||^_ 1 ||w; m ||^_ 2 +C , (a 2 )||F m ||^_ 1 ||u m ||2 /s 

Finally, the estimates for J, II, and III allow us to conclude 

"77 1 Pm 1 1 #3-2 < (— A + e)||u; m || HS -i + C||F m || HS _i||ti; m || //s _2 

+C ||'Um||H s - 2 |p m ||Hs-l||'U; m ||#s-2 + C ||-rm||#s-2||w TO ||#s-i||M m ||# J! 

-\-C I \F m \ \ H s-i I \u m \ \ H s + C\ \F m \ \jjs-i | \u m \ \jjs-i. (3.2.18) 

Using inequality (13.2.111) and the elliptic regularity estimates of Proposition 13.2. 1\ 
we can bound the right hand side of (13.2. 18f) by only t, \\F \\ H k, and ||w m || ff fc for 
k < s — 1. This is precisely the bound we need to show the existence of classical 
solutions using Proposition 13.2.11 As an example of this computation, we include the 
case when s = 4. 
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Example s = 4. We now use Proposition 13.2.11 and inequality (j3.2.1ip to write the 
estimate above for s = 4, in terms of w m and F alone. We begin with the second 
term, 

\\Fm\\ 8 H 3\\w m \\ 2 H 2 < C\\F \\ 8 Ha (l + 8t\\u m \\ H 4 + O(t 2 ))\\w m \\ 2 H2 

< C||F |||3(H-O(t 2 ))|W|^ +tC||Fo|^3||« m ||^IK||^ 2 

< C(t, \\Fo\\ H 3)\\w m \\ 2 H 2 + CdlFoWns^lwmW^ + C(t, \\F \\ H 3)\\u m \\ 2 H 4 

< C(t,\\F \\ H a)\\w m \\ 2 H 2 +C(t, I \F \ M\ \w m \ \%2, 



where we have used Young's inequality for the third inequality and inequality (13.2.61) 
for the last. The third term in inequality (13.2.181) can be estimated by a repeated use 
of Young's inequality, 



|'«m||/HKm||/H|^m||H 2 < C\\w m 

< C\\w m 

+C(t, 

< C\\w m 



2 xr'ii,/ II 4 



2 H 2 + C{t,\\F \\ H B)\\w m \\i 2 

I -Pol \h3)\ \w m \\%i\\w m \ \ 2 H 2 + C(\\F Q \\ H s)\\w m \\ 8 H i 

H2 + C(t,\\F \\ H 3)\\w m \\%2 

\Fo\\ H s)\\Wm\\ 8 H i. 



The fourth term of inequality (13.2.181) can be estimated by 

\\F m \ \h 2 \ \ U rn\\H s \ \ u m\ \h 4 — C(t)\ 1-^0 1 Iff 2 I l M m| l# 4 + ^1 l-^ol Iff 2 I \ U m\ \h 3 I l M m| \ff 4 

< C(t, I \F \ \hs) I \u m \\ 4 m + C(t, | \F \\ H 2)\\u m \\% 3 

< C(t, | \F \ \ H 2) | \w m \\ A Hl + C(t,||F|| H2 )||«; m ||| 1 
+C(t,\\F \\ H 3)\\w m \\ 4 H 2, 
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by inequalities 03.2.5p . (13.2.6J) . and Young's inequality. The second to last term of 
inequality (13.2.181) 

||-^m||jy3|Pm||H 4 — ^(^)ll-^o||_f/ 3 l \ u m\ \h 4 + £C1|F o\\h 3 \ \U m \ \ h a 

< C(t,\\F \\ H 3)\\u m \\ 2 H 4 + C(t,\\F \\ H 3)\\u m \\%4 

< C(t,\\F \\ H 3)\\w m \\ 4 Hl + C(t,\\F \\ H 3)\\w m \\ 8 Hl 
+C(t,\\F \\ H s)\\w m \\ 2 H 2 +C(t, \\Fo\\ H 3)\\w m \\%2, 

where we have used Young's inequality to write ||it m ||#4 < C|| u m||#4 + C|l M m||#4- 
The final term of inequality (13.2. 18j) can be estimated by inequality (I3.2.3P as 

iip ii 2 lUi ii 4 <r r\\j? ii 2 nj.r'iip ii 16 ^ 2 iu„ n 4 

<r r\\\T? M 2 j-IIp M 18 _i_iif M 34 liu„ M 4 

_i <-> |JI m M.ff 3 "T \\ r m\\H 3 ' \\ r m\\H 3 i 1 1 "'m 1 1 H 1 

< C(t, \\Fo\\ H 3)\\w m \\ 4 H i +tC(\\F \\ H 3)\\u m \\ H A\ w m\\H 1 

< C(t,\\F \\ H 3)\\w m \\* Hl + C(t,\\F \\ H 3)\\w m \\ 8 m 

+C(t, \\F \\ H 3)\\u m \\ 2 H 4 

< c(t,\\F \\ H 3)\\wj\jp. + c(t,\\F \\ Ha )\\wj\ B H1 

+C(t, \\Fo\\H 3 )\\Wm\\H2. 

Finally, combining the estimates above we achieve the a priori estimate 

-jT I km | lip < {-\ + £)\\w m \\ 2 H3 +C\\w m \\ 2 H 2 + C\\w m \\ 4 H 2 + C\\w m \\%i + C||lU m ||^i, 

(3.2.19) 
where C = C(t, a, \\F \\h 3 ) > 0. 
Convergence to a strong solution to equation (13.1.11) . 

From inequality (13.2.191) and equation (13.2.131) . we may conclude that for some 
T > 0, w m is bounded in L°°(0, T; H 2 iv ) n W 1 ' 00 ^, T; L 2 div ) n L 2 (0, T; #f iv ) uniformly 
in m, and hence by inequality (13.2.61) . 



u m is bounded in L°°(0, T; H^) n ^°°(0, T; H 2 div ) n L 2 (0, T; H 5 div 
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uniformly in m. By the weak compactness theorem, there exists a subsequence u nik 
such that for all T > 0, 

u mk - u in L°°(0, T; H 4 div ) R W^O, T; tf d 2 iv ) n L 2 (0, T; tf d 5 iv ). 

Since H /1 ' oo (0, T; H\ v ) C C°([0, T]; H% iv ), Arzela-Ascoli compactness criterion implies 



u mk ^umC°([0,T};H? 



divJ 



Furthermore by interpolation, u m e VP°'°°(0, T; H% iv ) n W^°°(0, T; Hj iv ) implies 

n mfc ^ M inC 5 ([0,T];F d 4 r /), 
and therefore by Sobolev's embedding theorem ( 12.2.91) . 



u mk ^umC°{[^T]-C t 



2 > 
divy 



for 5 taken sufficiently small. Thus all the terms on the right-hand side of equation 
(13.2. 13D converge strongly and d t u m converges weakly. To see that u G C°([0, T]; H^ v ), 
it suffices to show that | \u(t, •) | \h 4 is continuous on [0, T], but this follows from the 
inequalities (|3~23D and 03.2.191) . 

The uniqueness of classical solutions to (13.1.11) follows the same arguments made 
in [22] • The only deviation being the evolution of the difference of two solutions of 
(13.2.121) is given as 

d t (u 1 -U 2 ) + (Pf - P 2 a )Vl«! ■ Ml + ? 2 a V«l • («l - « 2 ) + ? 2 a V( Ml - M 2 ) ■ M 2 

= up? div((Fi - F 2 ) ■ V«i) + up? div(F 2 • Vui) + vP 2 div(F 2 ■ Vu 2 ) 
= uP? div((F a - F 2 ) • Vmi) + UP? div(F 2 • V(«i - u 2 )) 
+z/(Pf-P 2 a )div(F 2 -Vn 1 ), 

where F± solves equation (13.1. left with u = Ui for % = 1,2. That is, if 

y(t) := |K(£) - u 2 (t)\\m + \\Fi(t) - F 2 (t)\\%i 
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we use the fact that u E L°°(0, T; H^ lv ) fl L 2 (0, T; H^ 1 ) and computations similar to 
that used for existence to obtain the differential inequality 
d 



d -y(t) < C(t)y(t), C(t) = C{\\ui\\ H s, H^HhO, i = 



5 -"5 



from which uniqueness follows. I 

Remark 3.2.3 As a — > 0, the anisotropic LANS-a equations should reduce to the 
Navier-Stokes equations. Once it is shown that F — ► Id as a — ► 0, the anisotropic 
equations reduce to the isotropic equations. Marsden and Shkoller WT^ have proven for 
s > 3, solutions to the isotropic LANS-a equations converge in H s for short time on 
intervals which are governed by the existence theory for the Navier-Stokes equations. 
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Chapter 4 



Numerical solutions to the 



anisotropic LANS-a equations 



4.1 Introduction. 

Examining the behavior of fluid flow in elementary domains is more than a common 
numerical test for a new fluid model, it is also a means of highlighting interesting phe- 
nomena inherent in the model. The direct numerical simulation (DNS) of turbulent 
flow at small to moderate Reynolds number has become an important computational 
tool in understanding large scale turbulence motion. Unfortunately, DNS are still 
computationally expensive in turbulent regimes. Alternate approaches to brute-force 
DNS are the Reynolds averaged Navier-Stokes (RANS) simulations and large eddy 
simulations (LES). The RANS equations (I1.2.3P are briefly discussed in Chapter [TJ 
In LES, a spatial averaging operator (filter) is applied to the Navier-Stokes equations 
to obtain a new set of equations for the averaged (filtered) variables. Due to the 
computational limit of DNS for large Reynolds number, LES have become one of the 
standard methods in solving for fluid flow. 

The behavior of small spatial scales in turbulent flow is often characterized by 
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statistical isotropy and homogeneity away from the boundary of the fluid container. 
Therefore, the isotropic LANS-a equations appear to be an appropriate model in 
isotropic regimes and have been studied recently from the numerical point of view. 
As a proposed model for large scale turbulence, Mohseni et. al. [21] compare the 
isotropic LANS-a equations to known results for DNS and LES methods. The authors 
demonstrate the utility of the isotropic LANS-a as a sub-grid stress model for three- 
dimensional isotropic forced and decaying turbulence. They perform two sets of 
forced isotropic turbulence simulations and compare their results to DNS and LES 
results where appropriate. In the LANS-a simulations, Mohseni et. al. conclude 
that selecting an appropriate a is a delicate comprise between the accuracy of the 
model at the the large scales and the minimum resolution requirements. In particular, 
the LANS-a equations accurately mimic the behavior of the Navier-Stokes equations 
at large spatial scales as long as a minimum resolution is observed. The accuracy 
improves as a — > 0, but the computation requires higher resolution. The higher 
resolution comes with a price: an unresolved LANS-a computation could result in 
the contamination of large scales and a loss of accuracy. 

In 1998, Chen et. al. [5] studied the mean velocity of turbulent channel and 
pipe flows. They proposed using stationary solutions to isotropic LANS-a equations 
as a closure approximation for the Reynolds-averaged equations. Since in the near- 
wall region the fluctuations are highly anisotropic, their results were only in good 
agreement with experimental data away from the viscous boundary layer. 

By restricting the anisotropic LANS-a equations to the channel, Coutand and 
Shkoller [11] propose a turbulent channel theory that models the large scale fluid 
motion throughout the entire domain. Unlike the isotropic equations, the solutions 
to the anisotropic equations consist of both a mean velocity field u and a 3 x 3 co- 
variance tensor F, defined, in detail, in the next section. An important property of 
F is its degeneracy to zero at the wall. Coutand and Shkoller [11] show that near 
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the wall, the degeneracy scales like dy/\ \ogd\, where d is the normalized distance 
function to the wall. The authors compensate for the degeneracy in the bound- 
ary layer by working in weighted Sobolev spaces. In this functional framework, they 
prove the global-in-time existence and uniqueness of weak solutions to the anisotropic 
LANS-a equations. In particular, they restrict the anisotropic equations to the chan- 
nel and make the assumption that the initial covariance tensor is given in the form 
F = F(0,x,y,z) = p(z) Id. By assuming the fluid is moving in one direction, the 
anisotropic LANS-a equations reduce to a one-dimensional partial differential equa- 
tion for the mean velocity. The authors use the Galerkin method to obtain the 
existence of weak solutions to the anisotropic equations. 

In this chapter, we use the anisotropic LANS-a equations as a numerical model 
for two classical examples of laminar velocity profiles. Unlike other approaches to the 
modeling of tubulence, the Lagrangian averaging approach allows us to capture the 
large scale motion of a fluid in laminar regimes. The first example that we consider is 
the Poiseuille flow in the channel and pipe domains. Laminar Poiseuille flow occurs 
when an incompressible fluid with no-slip boundary conditions is driven by a constant 
upstream pressure gradient, yielding a symmetric parabolic stream-wise profile. We 
assume the velocity u satisfies the steady Navier-Stokes equations and we numerically 
solve the anisotropic equations for the initial covariance tensor F . In particular, we 
calculate the matrix F such that pair (u, F) is a solution to the anisotropic LANS-a 
equations. The degeneracy rate of our numerical solution F near the wall of the 
channel is in good agreement with the logarithmic decay rate given by Coutand and 
Shkoller |llj . In addition, we show that in the boundary layer, \\F(t, -)IU°°(n) does 
not remain bounded as t — > oo, answering, at least numerically, a question stated in 

nu. 

In the last section of the chapter, we study shear flow solutions to the anisotropic 
LANS-a equations. Shear flow occurs in the channel when one side of the boundary 
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is moving while the other side remains fixed. The derivation by Marsden and Shkoller 
[22] relied upon no-slip boundary conditions for the mean velocity field u. In order 
to find a solution to the inhomogeneous problem, we use the classical technique of 
introducing a new vector field that vanishes on the boundary and that solves the 
anisotropic equations with additional forcing terms. For the domains and the initial 
data that submit a unique solution to the new formulation, we are able to solve for the 
unknown mean velocity. In Section 14.3} we show that shear flow velocity solutions to 
the anisotropic LANS-a equations exist if the initial covariance tensor is not required 
to be positive. This turns out to be a natural assumption and we compute the shear 
flow solutions directly. 

4.1.1 The covariance tensor F(t, x). 

The matrix F is defined by Marsden and Shkoller [22] as the ensemble average of the 
tensor product of the Lagrangian fluctuation vector. Specifically, if r\ and rf are the 
particle trajectories of the fluid velocity u and averaged fluid velocity u e , respectively, 
then the Lagrangian fluctuation vector £ e := if o r/ _1 (see Figure H~TT) and 

d 

e=0 

where < • > denotes the average over all possible solutions u £ (see [22]). Because 
these fluctuations are necessarily zero along the boundary, 

F(t, x) = for t > 

on the boundary of the domain. As noted in [TTj . it is unknown a priori whether 
1 1-^(^5 ')lli°° remains bounded for all time. In next section, we show numerically that 
I \F{t, ')\\l°° is strictly increasing in time in the viscous boundary layer of the channel 
and pipe under a steady mean flow. 



de 



e=o de 



e 
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y = r)(t,x 




e 



— ff o r\ 



• f = A. I te 

• F = 0ondtt 

lug = u + ew 

a; = 77(0,3;) ^(0,x) = ^(0,x) 

Figure 4.1: The construction of i<\ 

4.2 Homogeneous boundary conditions. 

Laminar Poiseuille flow occurs when an incompressible fluid in a straight channel, 
or pipe, is driven by a constant upstream pressure gradient, yielding a symmetric 
parabolic stream-wise velocity profile. In this section, we begin with the steady 
Poiseuille flow associated with the Navier-Stokes equations. We construct an initial 
covariance tensor such that the steady flow solution and the covariance tensor solve 
the anisotropic equations in the channel. First, we reduce the anisotropic LANS-a 
equations to the channel and pipe under the assumption that the initial covariance 
tensor is a multiple of the identity matrix. 



4.2.1 Given the initial covariance matrix. 

Channel. The three-dimensional channel is given by Q = M 2 x [—h, h] with coordi- 
nates x = (x, y, z). We shall assume that the velocity vector is of the form 

u(t,x,y,z) = (u(t,z),0,0), 

and that the initial covariance matrix is written as 



F = F{0,x,y,z) = p{z)Id, 
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such that p(±/i) = 0. 

Suppose p(z) is given. As discussed in [IT], when no-slip boundary conditions are 
prescribed for the mean velocity, the nonlinear term vanishes and equation (13.1.11) 
reduces to the following system 

d t (u - a 2 (puj) - v((pu'y - a\p(puj)>) = -c, 

(4.2.1) 

u(t, ±h) = 0, u(0, z) — Uq(z), 
where c := d x p is constant since d y p = d z p = 0. To address the degeneracy at the 
boundary, Coutand and Shkoller define a weighted Sobolev to set as the functional 
framework. This space is based on the following function: For e < 5 << 1 define 
Pcs( z ) £ C°°[— h, h] to be the positive function 

f d^/\h^d\, < d{z) < e, 
Pcs{z) = < 

( 1, d(z)>8, (4.2.2) 

d = normalized distance function to the boundary. 

Coutand and Shkoller prove the existence of a unique global weak solution to equation 

(14.2.11) in a weighted Sobolev space. In addition, they show that the covariance tensor 

F must degenerate like pes i n the viscous boundary layer. In the next section, we 

verify that our numerical solution p matches pes near the walls of the channel. 

Pipe. The three-dimensional pipe is given by Q = M. 3 with coordinates x = (x, y, z) 

such that < \Jx 2 + y 2 < a. We shall assume that the velocity vector is radially 

symmetric of the form 

u(t,x,y,z) = (0,0,u(t,r)), (4.2.3) 



where r := y ' x 1 + y 2 and that the initial covariance matrix 

F = F(0,x,y,z) = p(r)ld. 

We let r\ denote the Lagrangian flow of u satisfying the initial value problem (11. 2. 5ft . 
Then definition (I4.2.3P implies that 

r/(t, x) = yx,y,z-\- I u(s, r) ds J . 
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Therefore 



F(t,x) := Dn(t,y) • F (x) ■ Dr)(t,x) T = p{r) 



i o m 

T 

1 y -u 



^U y -U l+U 2 

r r 



(4.2.4) 



where U(t, x) := / d r u(s, r) ds. Together with definition (14.2.31) . C\i : = div(Vu ■ F) 
Jo 



reduces to 



Cu = (0, 0, - (rpuj) 



with u'(t, r) := d r u(t, r). The nonlinear term in the anisotropic LANS-a vanishes and 
the equations reduce to the following system 



8 t (u - a 2 -(rpu'Y) - - ((rpu'Y - a\rp(-(rpu')')')' 



(4.2.5) 



u(t,a) = 0, u(0,r) = uo(r), 
where c := d z p is constant since d x p = d y p = 0. 



4.2.2 Given the mean velocity profile. 

In this section, we provide the numerical results for steady channel and pipe flow 
assuming no-slip boundary conditions. We start with the steady solution u to the 
Navier-Stokes equations and proceed to find a p, and therefore a F, such that (u, F) 
solve the anisotropic LANS-a equations. 

Channel. We begin by assuming the velocity field is the steady flow given as 
u(t, x,y, z) = (u(z), 0,0) that satisfies the no-slip boundary condition u(±h) = 0. 
For the Navier-Stokes equations, the classical Poiseuille flow is given by 



u(z) = (3(h 2 -z 2 ) 
p(x) = Po — 2p(3x, 



(4.2.6) 



where /3,po > 0. In order for (u, p) to solve equation (14.2. ip with the velocity u(z) 
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given by equation (14.2.61) . the function p(z) must solve the ordinary differential equa- 
tion 

(z P y - a\p(z P yy = i, 

p(±h) = 0. 
Solving equation (14.2.7J) numerically leads us to the following Proposition. 



(4.2.7) 



Proposition 4.2.1 Let the function u be the steady solution . 2. 6] ) to the incom- 
pressible Navier-Stokes equations in the smooth three-dimensional channel. The pair 
(u, p) solves the channel anisotropic LANS-a equation fi4.2.1\ ) when the function p(z) 



solves equation ^4.2.1 ). A numerical solution is given in Figure 4 -2(a) when a = 0.1 
and h = (3 = 1. In addition, the covariance tensor F is defined as 

1 + A(3 2 t 2 z 2 -20tz 
F(t J z):=p(z)D V (t,x)Dr ] (t,x) T = p(z) 1 

-2(3tz 1 

Remark 4.2.1 We assume that p(z) is positive and shares the same symmetry across 
the channel as the velocity u(z). In addition, we assign in the center of the channel 
the values p(0) = A and p'(0) = 0. By a numerical shooting method, we find a value 
for A that forces the zero boundary condition p(l) = 0. This solution is plotted in 



Figure 4- 2(a). 

Remark 4.2.2 The numerical solution p decays at the same rate as the function 



pes near the boundary. In Figure 4- 2(b), we plot both functions near the wall of the 
channel. 

To illustrate the dynamics of F in the channel, we calculate its eigenvalues and 
conclude the following Proposition. 

Proposition 4.2.2 The eigenvalues ofF(t,z) are 

Ai(M) = p{z), 

A 2j3 (i,z) = \p{z) (2 + U 2 ± \U\VW + t) = p{z) (l + 2(3h 2 z 2 ± 2(3t\z\ ^1 + (3H 2 z 2 
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The eigenvalues are plotted in Figure 4-2(c) for the case of a = 0.1, (3 = 1, and the 
function p given as the solution to equation ft4.2.7\ ). Furthermore, \\F(t, -)IU°°(n) = 
0(t 2 ) and although our velocity is steady throughout the channel, the covariance tensor 



\F(t 



i-)\\L°°(n) 



oo as t — > oo 



near the boundary in the direction of the flow. This is illustrated in Figure 4- 2(d). 
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2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 




0.9 0.91 0.92 



(a) The solution p(z) of equation 14. 2.71 . 



(b) The solution p(z) and the function p(z) used in 
1221 near the boundary. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 




(c) The eigenvalues of F(t, z) at t = 2. 



(d) The evolution of ||F||i= 



Figure 4.2: For h = /? = 1 and a = 0.1. Because of symmetry, the graphs (a), (c), and (d) arc on 
[0, 1] only. 
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Proof. The eigenvalues are straightforward to compute. Since the largest eigenvalue 
A 2 (-,x) = 0{t 2 ) and A 2 < ||F(*,*)||£«», then \\F(t, -)\\ L oo = 0{f). The vector 

tz + VI + t 2 z 2 
v(t,z) = C 

-1 

solves the eigenvalue equation Fv = X2V. Therefore, as t — » 00, the norm of F 
increases to infinity in the direction of the flow (x-direction). | 

Pipe. The format and conclusions of this subsection follow closely the results in the 
channel. We begin by assuming the velocity field is a steady velocity field given as 
u(t, x) = (0,0, u(r)) that satisfies the no-slip boundary condition w(±a) = 0. For the 
Navier-Stokes equations, the classical pipe flow is given by 

u(r) = (3(a 2 - r 2 ) 
p(r) = po — 4i//3r, 

where /3,po > 0. In order for (u, p) to solve equation (14.2.51) with the velocity u{r) 
given by equation (j4.2.8p . the function p{r) must solve the ordinary differential equa- 
tion 



(4.2.8) 



r(rp' + 2p) - a 2 (rp(rp" + 3p'))' = 2r 
p{a) = 0. 
As in the channel domains, we arrive at the following Proposition. 



(4.2.9) 



Proposition 4.2.3 Let the function u be the steady solution fiJJTB ) to the incom- 
pressible Navier-Stokes equations in the smooth three-dimensional pipe. The pair 
(u, p) solves the pipe anisotropic LANS-a equation fl^.£.<5p when the function p(r) 
solves equation (ZJTW . A numerical solution is given in Figure \J7^ when a = 0.1, 
(3 = 1, and the radius of the pipe is a = 1. In addition, the covariance tensor F is 
given by \J^.2.J^ and has the same eigenvalues as the channel covariance tensor. 
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0.7 0.8 



Figure 4.3: The solution p(r) of equation (J4.2.9I ) 



4.2.3 Conclusions. 

As expected, the anisotropy of the fluid is of fundamental importance in bounded 
domains. We have demonstrated that the dynamics of F are just as important in 
laminar regimes as they are in turbulent regimes. Supposing that the covariance 
tensor is a constant multiple of the identity matrix, while accurate for very short 
time and in the center of the channel, becomes an inaccurate assumption in the 
viscous boundary layer. In both the channel and the pipe geometry, the evolution 



of F, illustrated by Figure 4.2(c) , is non-decreasing in time and achieves its greatest 



value in this boundary layer. As time increases, the location of the maximum value 
of ||.F||l°° in the channel approaches the limit z = 0.89. Therefore an accurate model 
for fluid motion in the entire channel or pipe should be founded on the anisotropic 
model, rather than the isotropic version, in this boundary region. 

The logarithmic degeneracy rate of the covariance tensor in the channel agrees 
well with the decay rate of the function computed in [TT] . Since the dynamics of fluid 
motion in the channel and pipe are similar, it is not surprising that the covariance 
tensor for the channel and pipe are the same. The next step in studying the numerical 
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properties of the anisotropic model is to solve for the covariance tensor when the 
velocity is given as a time-dependent solution to the Navier-Stokes problem. This is 
a necessary step in understanding the dynamics of the covariance tensor F. 

4.3 Inhomogeneous boundary conditions. 

In many physical models, ranging from the study of blood flow to the modeling of 
earthquakes, at least one of the boundary components of the fluid container are in 
motion. In these cases, the mean velocity of the fluid at the wall will no longer be 
zero. In this section, we study the anisotropic LANS-a equations with inhomogeneous 
boundary data. As a specific example, we consider the mean fluid motion in a channel 
when one of the boundary walls is not fixed. This motion is called shear flow. 

4.3.1 The LANS-a equations with inhomogeneous boundary data. 

Let Q be a bounded fluid container in IR n with boundary dQ and suppose that the 
mean velocity field u = g ^ on dQ. To find a solution pair (u, F) that solves the 
anisotropic LANS-a equations with inhomogeneous boundary conditions, we choose 
a divergence-free vector field v (£, •) G C°°(Q) such that v = g on the boundary. With 
v given, we define a new vector field w = u — v. The vector field w is zero on the 
boundary and solves the anisotropic LANS-a equations with u replaced with w + v. 
Namely, we now search for a solution to the following system of partial differential 
equations 

(1 - o?C) (d t w - uFCw + div(> <g> w + w ® v)) = - gradp - (1 - a 2 C) (d t v - vFCv) 

div w(t, x) = 0, w = on dQ, 

d t F + VF-w- (Vw ■ F + [Vw ■ F} T ) =VF-v + Vv-F+ [Vv ■ F} T , 

w(0, x) = u (x) — v(0, x), F(0, x) = F (x). 

(4.3.10) 
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In the three-dimensional torus T 3 , the trivial existence and uniqueness of a solution 
to (14.3. 10p follows from the smoothness of v, Theorem 1 in [22], and Theorem 13.2.11 
in Chapter [31 This result is stated as the following proposition. 

Proposition 4.3.1 For s > 7/2, and w G #j iv (T 3 ), F G [#* cr (T 3 )] 3x3 ; with F > 
0, there exists a unique solution (w, F) with w G C°([0, T]; H^ iv ) fl L 2 (0, T; H^ 1 ) and 
F G C°([0,T]; [Ft*] 3 * 3 ) to equations \J^.3. 1U\ ), where T depends on the initial data. 



In arbitrary bounded domains, it is unknown whether solutions exist to the 
anisotropic equations. As we demonstrate in the next section, under certain lim- 
iting conditions we can find a solution to equation (14.3. 10p . 

4.3.2 Given the initial covariance matrix. 

To study shear flow velocity solutions to the anisotropic LANS-a equations, we need 
to restrict the full equations to the three-dimensional channel by making a num- 
ber of limiting assumptions. We assume that the initial covariance tensor Fq = 
F(0,x,y, z) = p(z)Id is given and the mean velocity has the form u = (u(t, z), 0, 0) 
for z G [—h, h], where u^Oon the boundary. 

As was done in the general case above, we choose a vector field v(t, •) G C°°[— h, h] 
such that v (t, ±h) = u(t, ±h) for all t > 0. The vector field w = u — v then solves 
the following one-dimensional problem 



(4.3.11) 



C a w = -c + f, 
w(±h) = 0, w(0, z) = Uq{z) — v(Q, z) 

where C a is the linear operator defined as 

C a u := d t (u - a 2 (pu')') - v((jmj - a 2 \p{pu ')")'), 

and / := —C a v. Then u(t,x,y,z) := (w(t,z) + v(t, z),0,0) solves equation (14.2.11) 
with inhomogeneous boundary conditions. 
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4.3.3 Given the mean velocity field. 

Suppose that we are given the steady velocity vector field u = (u(z), 0, 0) with u(z) ^ 
on the boundary. To find a solution to the anisotropic model, we need to find a 
function p(z) such that 

u((pu f y - a?{p{pu')")') = -c. 

But since the vector field u doesn't satisfy the no-slip boundary conditions, we do not 
know the correct boundary conditions for p. Rather, we choose v G C°°[— h, h] with 
v = u on boundary Then w = u — v and the pair (w, p) solve 

C a w = -c + f, (4.3.12) 

where / := —C a v. Since w = on the boundary we also have the boundary condition 
p(±/i) = 0. Since w is given we may find p satisfying equation (14.3. 12p with zero 
boundary conditions. However, unlike the case when u satisfies no-slip boundary 
conditions, the forcing in equation ( 14.3.121) removes any a priori statement about the 
positivity of p. 

Shear flow solution. As an example, we show the existence of a shear flow velocity 
solution to the anisotropic LANS-a equations. The steady shear flow velocity solution 
to the Navier-Stokes equations with the boundary conditions u(—h) = and u(h) = 
U >0 is 

u(z) = ^-(z + h) (4.3.13) 

with the pressure function p = 0. The function 

w{z) = -^(h-z)(h + z) 

is zero on the boundary and solves (14.3.121) with / = — C a I -rr^{z + h) 2 J when p(z) 
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is a solution of 

p'-a\pp")' = 0, 
p{±h) = 0. 



Namely, 



2 



±p(z - a)(z -b) for -h<a< z<b<h 



otherwise 



If we want p £ C 1 [— /i, /i] then a = —h, b = h. We conclude that the pair (u, p) with 

1 
2c7 2 
which solves equation ( 14.2.11) with inhomogeneous boundary conditions. 



u define by equation (14.3.131) and p(z) = — -(z — h)(z + h) is a shear flow solution 
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